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1. 


INTRODUCTION 


This  final  report  describes  research  into  high-altitude  infrared 
(IR)  radiance  models  done  under  Contract  No.  F19628-87-C-0130.  A  major 
result  of  this  effort  has  been  the  development  of  predictive  models 
which  have  been  combined  into  a  general  calculational  model  for 
atmospheric  radiance.  These  models  include  sources  of  high-altitude 
NLTE  (Non  Local  Thermodynamic  Equilibrium)  radiation  in  the  IR  spectral 
region.  Many  of  the  models  have  been  combined  into  SHARC,  the  Strategic 
High  Altitude  Radiance  Code.  SHARC  calculates  atmospheric  radiance  in 
the  2-40  fjtn  spectral  region  for  arbitrary  line-of-sight  (LOS)  paths  that 
lie  above  50  km,  including  space  viewing.  The  latest  version, 
SHARC-2,^'2  also  includes  a  fully  integrated  auroral  model. 

At  lower  altitudes  molecular  collisions  are  many  and  frequent,  so 
collisional  processes  dominate  radiative  ones,  and  molecular  excited 
states  are  quenched  before  they  can  radiate.  The  molecules  are  then 
said  to  be  in  Local  Thermodynamic  Equilibrium  (LTE),  and  the  local 
temperature,  total  number  densities,  and  molecular  mixing  ratios  suffice 
to  determine  excited  state  populations.  At  higher  altitudes  (above 
-40-60  km) ,  competing  processes  like  radiative  decay  become  more 
important  and  eventually  dominate  collisional  deexcitation. 2 ^ 3 
Atmospheric  radiation  then  arises  from  emissions  by  molecules  that  are 
in  NLTE.  The  work  described  in  this  report  and  the  resultant  computer 
models  address  the  NLTE  regime. 

The  NLTE  models  and  the  resultant  SHARC  Code  are  described  in 
Sections  2  and  3.  In  related  activities,  data  taken  under  the  SPIRIT, 
EXCEDE,  and  SPIRE  programs  have  been  analyzed  and  compared  to  model 
predictions.  The  SPIRIT  1  data  were  used  to  study  O3  emission  in  the 
9.6  pm  band,  determination  of  O  and  H  atom  profiles,  and  NO  and  CO2 
emissions  in  their  5.3  and  15  pm  bands.  The  EXCEDE  flight  in  April, 
1990,  returned  much  useful  data;  Section  5  presents  an  initial  look  at 
the  spectral  data.  The  SPIRE  data  were  used  to  help  validate  parts  of 
SHARC.  Other  modeling  and  concepts  for  possible  new  measurements 
approaches  are  discussed  in  Sections  6,  7,  and  8. 
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THE  SHARC  COMPUTER  MODEL 


Under  NLTE  conditions,  both  the  ground-  and  excited-state  densities 
must  be  modeled  for  each  molecule  and  transition.  The  population  of  the 
various  excited  vibrational  states  deviate  most  strongly  from  LTE.  The 
relevant  NLTE  processes  fall  into  three  groups:  radiative,  collisional, 
and  chemical: 


•  Radiative:  Pumping  mechanisms  include  absorption  of 

sunshine,  upwelling  earthshine,  and  photons  emitted  by 
other  molecules.  Subsequent  spontaneous  emission  then 
results  in  the  observed  radiation. 

•  Collisional:  Non-reactive,  inelastic  collisions,  which 

result  in  vibrational  and  rotational  energy  exchange  or 
transfer,  are  still  significant  in  the  upper  atmosphere. 
Although  collisions  are  less  frequent  due  to  low  number 
densities,  molecule-radical  interactions,  such  as 

interactions  with  O  atoms,  and  near-resonant  interactions, 
such  as  V-V  exchange  between  N2  and  CO2,  are  important. 

•  Chemical ;  Inelastic  two-  and  three-body  collisions 

combined  with  chemical  reactions  frequently  leave  highly 
excited  products  which  subsequently  cascade  down  to  the 
ground  state  by  emitting  photons.  An  example  is  ozone 
recombination^  via  O  +  O2  +  M  -*  O3*  +  M  and  formation  of 
vibrationally  excited  NO  via  O2  +  N(2d)^  (atomic  nitrogen 
is  formed  during  auroral  activity). 


Identifying  and  quantifying  the  important  IR  radiators  depend  on 
data  obtained  from  field  measurements;  these  indicate  that  the  major 
sources  (and  bands)  above  80  km  are  NO  (5.3  pm),  CO2  (15  pm  and  4.3  pm), 
and  O3  (9.6  pm).  Other  emitters  of  varying  importance  include  CO 
(4.6  pm),  O2  (1.27  pm),  and  OH  (3.0  &  1.5  pm).^'^  At  lower  altitudes, 
H2O  (2.7  &  6.3  pm)  and  OH  become  increasingly  important,  and  the 
strength  of  NO  emissions  decreases. 

Various  models  that  describe  these  processes  have  been  synthesized 
into  an  overall  computer  model,  SHARC  (the  Strategic  High-Altitude 
Radiance  Code). SHARC  calculates  NLTE  and  LTE  radiance  and 
transmittance  for  paths  between  50  and  300  km  (including  space  viewing) 
within  the  2-40  pm  spectral  region  and  with  a  resolution  of  0.5  cm"^. 
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It  allows  arbitrary  paths  (e.g.,  limb,  horizontal,  vertical,  or  slant) 
between  any  two  points  within  or  without  the  atmosphere.  The  code 
presently  includes  the  six  most  important  ambient  molecules  at  high 
altitude:  CO2,  O3,  NO,  H2O,  CO,  and  OH.  In  addition,  SHARC  calculates 
auroral  emission  from  CO2,  NO,  and  NO"*". 


SHARC 

is 

a  completely 

new 

code. 

Modular  construction  has  been 

emphasized 

so 

that  models 

and 

model 

parameters 

can  be  modified 

or 

upgraded  as 

additional  data 

and/or  better  models 

become  available. 

It 

is  a  successor  to  the  AFGL  HAIRM  Code.  Both  SHARC  and  HAIRM  are 

based  on  analyses  of  various  AFGL  lab  and  field  measurements  for  these 
radiative  processes.  The  molecular  line  parameters  are  taken  from  the 
1986  HITRAN  line  atlasl^  and  supplemented  with  additional  lines  for  the 
high  vibrational  states  of  03(1/3)  in  the  10-12  pm  spectral  region^^  and 
lines  for  auroral  emission  by  NO  and  NO"*".^^  Parameters  in  the  SHARC 
data  base  can  be  upgraded  as  additional  data  are  analyzed.  This 
upgrade  ability,  its  integrated  auroral  model,  flexible  geometry 
package,  and  interactive  input  module  represent  significant  advances 
over  HAIRM. 

The  three  major  SHARC  components  are  indicated  in  Fig.  1.  The 
INTERPRETER,  which  is  a  modified  version  of  that  furnished  with  the 
Sandia  CHEMKIN  code,^^  prepares  the  differential  equations  for  the 
chemical/kinetic  reaction  scheme.  This  module  is  run  only  when  changes 
have  been  made  to  either  the  reaction  scheme  or  any  associated  rate 
constants.  The  SHARC  component  calculates  the  required  atmospheric 
transmittance  and  radiance.  The  chemistry  modules,  ambient  and  auroral, 
utilize  differential  equations  set  up  by  the  INTERPRETER.  The  PLOTTING 
package  generates  spectral  radiance  plots. 

2.1  The  INTERPRETER  Program 

Within  SHARC,  the  CHEMKIN  module  computes  the  steady-state  and 
time-dependent  number  densities  of  vibrationally  excited  atmospheric 
species  from  a  set  of  chemical  kinetics/reaction  mechanisms.  The 
chemical  kinetics  scheme  describes  molecular  formation,  all  forms  of 
vibrational  energy  transfer,  and  the  absorption  of  solar  and/or 
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Figure  1.  Major  Components  of  SHARC. 


earthshine  radiation.  The  INTERPRETER  program,  which  is  based  entirely 
on  and  includes  subroutines  directly  from  the  Sandia  Livermore 
INTERPRETER  code,!^  reads  a  symbolic  description  of  the  chemical 
kinetics  mechanisms  and  creates  a  binary  "linking”  file  with  all  the 
kinetics  information.  The  Sandia  CHEMKIN  package  is  described  as^5  -a 
general-purpose,  problem-independent,  transportable,  FORTRAN  chemical 
kinetics  code.”  The  SHARC  INTERPRETER  is  a  modified  version  from  which 
information  on  elements  in  the  periodic  table  and  the  thermodynamic  data 
base  (useful  for  combustion  reactions),  extraneous  to  the  atmospheric 
application,  have  been  removed.  The  INTERPRETER  reads  a  symbolic 
description  of  an  arbitrary  chemical  kinetics  mechanism  in  a  manner  that 
is  just  as  would  be  written  by  a  chemical  kineticist.  To  be  more 
specific,  consider  irreversible  kinetic  (energy  transfer  or  reactive) 
processes,  each  given  in  the  general  form 


K  ki  K 
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where  the  stoichiometric  coefficients  are  integers;  Cj^  is  the 

chemical  symbol  for  the  species;  and  kj^  is  the  rate  constant  for  the 

ith  process.  The  INTERPRETER  reads  these  equations  and  translates  the 
mechanism  into  the  appropriate  differential  equations 
I  K 

^  (**'ki~*^ki )  ^i  I  I  •  (2) 

i=l  k=l 

All  variables  defined  in  Eqs.  (1)  and  (2)  are  written  into  a  binary 
"linking"  file  for  use  by  SHARC.  Although  this  example  is  for 
irreversible  processes,  SHARC-2  allows  reversible  vibrational  to 
translational  processes.  The  INTERPRETER  only  has  to  be  run  once  for  a 
given  kinetics  mechanism  and  data  base. 

2.2  SHARC 

The  schematic  shown  in  Fig.  2  illustrates  the  calculational 

sequence  for  SHARC.  All  modules  are  called  by  the  MAIN  program.  The 

input  module  is  interactive  and  queries  the  user  for  parameters  needed 
by  the  rest  of  the  code.  The  chemistry  module  includes  two  components, 
the  chemical  kinetics  (CHEMKIN)  and  radiative  transfer  (NEMESIS) 
modules.  They  calculate  an  altitude  profile  of  populations  for  the 
various  molecular  excited  and  ground  states.  The  auroral  module 
calculates  auroral  enhancement  of  the  CO2,  NO,  and  NO'*’  ambient 

populations.  GEOMETRY  determines  the  line-of-sight  (LOS)  trajectory 
parameters  for  the  viewing  geometry  requested  by  the  user.  The  SPCRAD 
module  calculates  the  spectral  transmittance  and  radiance  for  the 
required  LOS  path(s). 

2 . 3  Ambient  Molecular  Populations 

The  first  component  of  the  SHARC  Ambient  module  is  CHEMKIN.  It 
computes  the  steady-state  number  densities  of  vibrational ly  excited 
atmospheric  species  for  each  atmospheric  layer.  Again  these  subroutines 
are  baaed  on  those  in  the  Sandia  Livermore  CHEMKIN  package. 
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Figure  2.  Calculational  Sequence  for  SHARC. 

CHEMKIN  uses  information  from  the  binary  "linking"  file  (species 
identification,  stoichiometric  coefficients,  and  rate  constant 
parameters)  to  construct  steady  state  equations  for  the  vibrationally 
excited  states 

I  K 

2<*'ki“‘^ki)'«i  J][  =  0  .  (J) 

i=l  k=l 


The  variable  definitions  are  identical  to  those  defined  in  Eq.  (2).  In 
general,  solving  Eq.  (3)  for  the  vibrationally  excited-state  densities 
requires  a  nonlinear  equation  algorithm.  However,  using  the  assumption 
that  the  of  the  vibrationally  excited  species  (i.e.,  the  unknowns  in 
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the  problem)  are  unity,  Eq.  (3)  simply  reduces  to  a  set  of  linearly 
independent  algebraic  equations  for  the  N  unknown  excited-state  number 
densities.  This  assumption  is  not  overly  restrictive  and  has  been  used 
previously  in  atmospheric  radiation  models. The  information 
necessary  to  solve  Eq.  (3)  is  simply  the  atmospheric  species  number 
densities  and  the  kinetic  temperature  from  which  the  rate  constants  are 
computed.  SHAHC  then  uses  a  lower/upper  (LU)  decomposition  algorithm  to 
solve  for  the  number  densities. 

The  second  component  is  NEMESIS,  the  Non-Equilibrium  Molecular 
Emission  and  Scattering  Intensity  Subroutine . It  computes  the 
enhancement  of  the  atmospheric  excited-state  layer  populations  due  to 
layer  radiative  self-trapping  and  layer-layer  radiative  pumping.  Some 
molecular  bands,  in  particular  the  4.3  pm  CO2  band,  are  optically  opaque 
to  emitted  radiation.  Photons  emitted  in  these  bands  may  be  absorbed 
and  emitted  many  times  before  either  escaping  the  atmosphere  or  becoming 
collisionally  quenched. 

Calculating  radiatively  enhanced  excited-state  populations 
involves: 

•  determination  of  layer  steady-state  source  populations, 
which  includes  excitation  by  external  light  sources,  sun 
and  earthshine,  and  molecular  collisions  and  de-excitation 
by  radiative  decay  and  collisional  quenching; 

•  determination  of  first-order  population  enhancements  using 
a  Monte  Carlo  simulation  of  source  photon  emissions  and 
their  subsequent  absorption  or  escape;  and 

•  determination  of  the  total  enhanced  populations  using  a 
recursive  orders-of-enhancement  approximation  which  is 
initialized  by  the  Monte  Carlo  first-order  results. 

The  source  populations  include  all  sources  of  excitation  exclusive  of 
the  internal  atmospheric  radiative  effects  of  layer  self-trapping  and 
layer-layer  pumping.  CHEMKIN  generates  the  source  populations  for  use 
by  NEMESIS.  The  Monte  Carlo  simulation  calculates  the  first-order 
enhancements  and  the  probabilities  that  a  photon  emitted  from  a  layer 
"i”  will  create  a  new  excited-state  in  each  layer  "j".  This  simulation 
involves  sampling  over  the  initial  emission  position,  emission 
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direction,  emission  frequency,  emission  line  strength,  and  length  of 
travel.  The  contribution  of  each  succeeding  order-of-enhancement  is 
determined  recursively  by 

(Pjt)  =  (Pk-l)(M)  .  (4) 

where  (P]^)  is  the  kth-order  layer  population  enhancement  matrix,  and  (W) 
is  the  layer-layer  absorption  probability  matrix.  Both  (P^)  and  (W)  are 
calculated  by  the  first-order  Monte  Carlo  calculation. 

The  end  result  of  each  NEMESIS  calculation  is  the  total 
excited-state  population  distribution  for  a  pair  of  vibrational  levels, 
where  the  lower  level  of  each  pair  may  itself  be  an  excited  level.  A 
cycling  procedure  between  CHEMKIN  and  NEMESIS  is  used  to  step  up  the 
vibrational  energy  ladder  until  the  top  of  the  ladder  is  reached. 
Finally,  these  enhanced  populations  are  passed  to  the  SPCRAD  module  for 
the  LOS  radiance  computation. 

2 . 4  Auroral  Molecular  Populations 

The  AURORA  module  performs  a  time-dependent  calculation  to 
determine  auroral  species  number  densities  as  a  function  of  altitude. 
The  Gear  method  for  stiff  differential  equations  is  used  to  numerically 
integrate  the  chemical  kinetic  rate  equations  that  describe  the  various 
auroral  chemical  and  energy  transfer  processes.  This  integration 
package  is  part  of  the  Sandia  CHEMKIN  codecs  thus  has  been 
extensively  applied  to  many  different  types  of  scientific  problems. 

Auroral  excitations  result  from  the  atmosphere's  stopping  of  high- 
energy  electrons  impinging  upon  it.  The  deposition  model  for  these 
primary  electrons  is  based  upon  the  work  of  Grun,  Rees,  and  Strickland 
and  has  been  discussed  in  some  detail  in  the  AARC  manual. 14  The  SHARC 
routines  for  energy  deposition  are  based  upon  those  used  in  AARC.  The 
ion  pair  production  rate  is  obtained  from  the  energy  deposition  rate  by 
assuming  that  35  eV  are  required  to  produce  each  ion  pair. 

Chemical  kinetics  mechanisms  which  lead  to  population  enhancements 
of  NO,  NO'*',  and  C02{4.3  pm)  are  included  in  the  SHARC-2  data  base.  The 
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relevant  chemical  reactions  and  energy  transfer  processes  starting  with 
the  initial  electron  deposition  are  also  described  in  the  AARC  rnanuall** 
and  in  Section  3.  The  mechanisms  include  kinetic  processes  for 
ionization,  electronic  and  vibrational  excitation,  recombination, 
chemical  reactions,  and  formation  of  the  secondary  electron  energy 
distribution. 

2 . 5  LOS  Spectral  Radiance 

Calculation  of  the  spectral  transmittance  and  radiance  for  the 
requested  LOS  path  uses  a  finite-difference  form  of  the  radiative 
transfer  equation.  A  single-line  equivalent  width  approximation  based 
on  the  Voigt  lineshape  is  used  to  calculate  segment  transmittances  and 
radiances.  This  approach  enables  the  spectral  radiance  to  be  calculated 
on  a  line-by-line  (LBL)  basis  at  a  spectral  resolution  as  high  as 
0.1  cm~l,  although  the  resolution  in  the  current  version  of  SHARC  is  set 
at  0.5  cm”l.  Since  the  radiation  computation  is  explicitly  based  on  the 
difference  of  the  upper  and  lower  state  populations,  it  is  equally  valid 
in  both  the  NLTE  and  LTE  regimes.  Line  strengths  and  locations  are 
taken  from  a  modified  line  file  generated  from  the  1986  AFGL  HITRAN  line 
atlas. The  major  modifications  were  addition  of  some  additional  lines 
and  changing  the  line  parameter  data  base  to  decompose  the  energy  of  the 
lower  state  into  vibrational  and  rotational  energies.  This  enables  the 
line  strengths  to  be  scaled  separately  for  the  rotational  temperature 
and  the  NLTE  vibrational  population  of  the  lower  level. 

There  is  a  trade-off  between  speed  and  accuracy  when  the  atmosphere 
is  layered  into  many  homogeneous  layers.  Times  for  LOS  radiance 
calculations  depend  linearly  on  the  number  of  layers.  Currently  the 
atmospheric  profiles  used  by  SHARC  are  layered  by  2  km  from  50  to  150  km 
and  by  10  km  steps  from  150  to  300  km.  Thus  the  LOS  path  for  a  limb 
calculation  with  a  tangent  height  of  50  km  traverses  129  atmospheric 
layers. 
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Initial  SHARC  validation  was  accomplished  by  comparison  to  field 
data  from  the  GL  SPIRE, SPIRIT  1,^®  and  FWI^®  programs.  In  addition, 
extensive  comparisons  have  been  made  to  other  GL  codes  to  verify  both 
the  performance  of  individual  modules  and  the  overall  code. 

The  SPIRE  measurements  were  taken  from  a  sounding  rocket  with 
spectrometers  that  had  a  spectral  resolution  of  about  2  per  cent.^^  The 
data  are  for  partial  to  full  limb  observations.  All  presented  data  have 
been  corrected  to  full  limb  radiances.  Figures  3  and  4  present 
comparisons  for  CO2  radiances  integrated  over  the  CO2  >^2  *^3  i>snda  at 
15  and  4.3  pm,  respectively.  The  two  calculations  in  Fig.  3  illustrate 
the  contribution  of  molecular  radiative  exchange.  The  NO-NEMESIS 
calculation  for  which  the  radiance  is  calculated  with  just  chemical  and 
earthshine  excitation  processes,  illustrates  the  importance  of  radiative 
excitation.  The  shoulder  around  100  km  arises  from  a  complicated 
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Figure  3.  Comparison  of  Calculated  and  Observed  Radiation  in  the 
15  pm  Band  of  CO2  for  Limb  Viewing  from  Space.  The 
NO-NEMESIS  Calculation  Does  Not  Include  Molecular 
Radiative  Excitation. 
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Figure  4.  Comparison  of  Calculated  and  Observed  Daytime  Radiation 
in  the  4.3  pm  CO2  Band  for  Limb  Viewing  from  Space. 


interplay  between  optical  opacity  of  the  primary  CO2  isotope  and 
excitation  of  C02(t'2)  ®  atoms.  Sunlit  SPIRE  scans  and  model 

calculations  are  shown  for  the  4.3  pm  t/3  band  in  Fig.  4.  Different 
phenomena  dominate  at  various  altitude  regimes.  Above  110  km,  direct 
solar  excitation  of  the  primary  isotope  dominates.  Below  this  altitude, 
the  radiance  is  strongly  enhanced  both  by  emission  from  the  minor 
isotopes  and  by  solar  pumping  of  the  2.7  pm  band  with  subsequent 
radiative  decay  of  the  hot  bands  in  the  4.3  pm  region.  The  strongest 
excitations  coupled  to  the  2.7  pm  band  are  the  and  1/2=2  states 

combined  with  additional  1^3  quanta.  The  shoulder  below  85  km  arises 
from  collisional  quenching.  As  the  density  increases,  these  NLTE 
excitations  are  quickly  quenched,  and  the  single  i/j  excitation  of  the 
primary  isotope  again  dominates  as  LTE  emissions.  There  is  a  small 
nighttime  shoulder  around  85  km  due  to  the  isotopes  (not  shown). 

Figures  5  and  6  show  data-model  comparisons  for  the  9.6  pm  O3  and 
5.3  pm  NO  bands,  respectively.  Both  model  calculation  were  done  for  a 
60°  North  (Subarctic)  Winter  Profile.  The  primary  source  of  the 
nighttime  limb  emission  for  O3  in  Fig.  5  results  from  three-body 
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Figure  5.  Comparison  of  Calculated  and  Observed  Nighttime 
Radiation  in  the  9.6  fjm  O3  Band  for  Limb  Viewing  from 
Space. 


e 


0) 

■o 

D 


CE 


Figure  6. 
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Comparison  of  Calculated  and  Observed  Radiation  in  the 
5.3  fjm  NO  Band  for  Limb  Viewing  from  Space. 
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recombination  of  O2  with  atomic  oxygen.  The  data  in  the  figure  are  from 
the  SPIRE  and  from  SPIRIT  1^®  field  measurements.  The  model  and  data 
are  in  good  agreement.  A  similar  comparison  to  emissions  in  the  NO  band 
is  shown  in  Fig.  6.  Agreement  with  the  SPIRE  data  is  quite  good,  while 
the  SPIRIT  measurements  are  about  a  factor  of  five  lower.  This  is 
indicative  of  the  effects  of  atmospheric  variability.  The  strength  of 
NO  radiation  is  proportional  to 

(OJ  (NO)  e"2600/T  ^  ^5^ 

where  T  is  the  atmospheric  temperature  and  [O]  and  [NO]  are  the 
concentrations;  thus,  emissions  in  this  band  have  a  strong  dependence  on 
the  ambient  number  densities  and  temperature.  As  discussed  in  Section 
4,  ambient  temperatures  during  the  SPIRIT  measurements  were  consistently 
lower  than  those  during  SPIRE. 


2 . 7  Multiple  Atmospheric  Regions 

Much  of  the  nomenclature  in  the  GEOMETRY  module  has  been  borrowed 
from  the  LOWTRAN-6  geometry  routines. 20  fj©  knowledge  of  LOWTRAN  6  is 
required  as  GEOMETRY  is  self-contained.  A  curved-earth  geometry  is  used 
for  all  LOS  path  calculations;  since  atmospheric  refraction  is 
insignificant  under  these  rarefied  conditions,  it  is  not  calculated. 

The  auroral  excitation  is  confined  to  a  finite  region  within  the 
atmosphere.  This  has  led  to  the  development  of  local  regions  with  their 
own  atmospheric  profiles.  A  local  region  can  receive  auroral  excitation 
due  to  a  prior  history  of  electron  dosing  or  be  just  a  different  ambient 
region.  Up  to  four  local  regions  are  allowed. 

Local  regions  are  defined  by  four  (latitude,  longitude)  points  plus 
the  upper  and  lower  altitudes  of  the  region.  The  four  (Iat,lon)  points 
must  form  a  convex  quadrilateral  which  has  no  interior  angles  greater 
than  180°.  Two  local  regions  must  not  intersect  one  another,  but  they 
can  share  common  boundaries.  If  a  pole  is  to  be  within  a  local  region, 
the  input  geometry  should  use  the  other  coordinate  system,  i.e.,  use 
polar  coordinates  when  the  local  region  includes  the  geomagnetic  pole. 
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Figure  7.  Layering  Schematic  with  Local  Auroral  Region. 

Usually  the  LOS  of  the  radiance  calculation  only  intersects  a  local 
region  twice.  A  typical  path  is  illustrated  in  Fig.  7.  However  four 
intersections  are  possible  for  a  limb  path  whose  minimum  altitude  is 
just  slightly  less  than  the  lower  region  boundary.  This  arises  because 
of  earth  curvature  effects.  For  excunple,  a  limb  path  at  89  km  combined 
with  a  lower  region  boundary  of  90  km  will  have  four  intersections  when 
the  horizontal  extent  of  the  region  is  several  hundred  kilometers  to 
each  side  of  a  long  tangent  path. 

2 . 8  Atmospheric  Spatial  Structures 

Natural  backgrounds  are  non-uniform  in  both  space  and  time  and  are 
subject  to  changes  due  to  geophysical  conditions  such  as  latitude, 
magnetic  activity,  and  localized  dynamic  effects.  Other  spatio-temporal 
non-uniformities  can  be  attributed  to  factors  like  aurora,  eddy 
diffusion,  winds  or  large  scale  circulations,  tides,  solar  wind,  and 
gravity  waves.  Structured  features  in  the  upper  atmosphere  can  be 
divided  into  three  general  classes: 

•  long-term  variations  of  the  quiescent  atmosphere, 

•  short-term,  wave-like  fluid  responses  to  forces  acting 
on  the  atmosphere,  and 

•  variations  due  to  energy  deposition  from  the 
magnetospheric  particles  and  fields. 
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Short-term  variations  are  the  primary  sources  of  "clutter"  for  systems 
applications,  that  is,  spatial  and/or  temporal  fluctuations  in 
atmospheric  radiance.  Long-term  variations  are  global  effects  which 
describe  average  properties  of  the  atmosphere  over  distances  larger  than 
those  related  to  an  observer's  FOV  or  LOS.  One  major  manifestation  of 
energy  deposition  is,  of  course,  auroral  activity.  Scales  associated 
with  some  typical  atmospheric  structures  are  given  in  Table  1. 


Table  1.  Atmospheric  Structure  Scales. 


•  LONG-TERM  VARIATIONS 

horizontal  scale 
vertical  scale 
temporal  scale 
occurrence 

•  SHORT-TERM  TURBULENCE 

horizontal  scale 
vertical  scale 
temporal  scale 

fractional  perturbation 

form  of  PSD 
occurrence 

•  SHORT-TERM  GRAVITY  WAVES 

horizontal  scale 
vertical  scale 
temporal  scale 
fractional  perturbation 
occurrence 

•  THERMAL  TIDES 

vertical  scale 
temporal  scale 
fractional  perturbation 
occurrence 


-  100  km 

-  scale  height  H  *  6  km  i 

-  1  hr 
seasonal, 

solar  activity,  etc. 


-  50  km 

-  12  hr 

^  =  20% 

P 

regularly,  due  to  solar 
heating  of  atmosphere 
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There  is  a  severe  lack  of  good  data  on  atmospheric  spatial 
structures  in  the  SO-300  km  regime.  Several  approaches  for  using  the 
basic  science  contained  in  the  present  SHARC  chemistry  and  radiation 
transport  modules  as  the  basis  for  a  systematic  phenomenological 
approach  to  clutter  were  studied.  Structures  with  large-scale 
variations  could  be  adequately  handled  for  many  applications  by 
expanding  the  number  of  default  atmospheres  to  include  latitudinal, 
diurnal,  seasonal,  geomagnetic,  and  solar-activity  effects.  These 
profiles  could  be  linked  to  an  available  model  like  MSIS,^^  that  has 
been  upgraded  to  include  profiles  for  important  IR  molecules  like  CO2, 
NO,  and  O3.  Small-scale  fluctuations  could  be  modeled  with  a  causal 
approach  in  which  spatial  structures  in  radiance  are  calculated  from 
spatial  variations  in  basic  atmospheric  parameters  like  density, 
composition  and  temperature.  These  structures  can  be  introduced  into 
SHARC  by  defining  different  local  atmospheric  regions,  each  having  its 
own  density  and  temperature  profile. 

This  approach  to  small-scale  clutter,  which  we  call  geometric, 
assumes  that  radiance  structures  result  from  atmospheric  regions  with 
different  temperatures  and  species  densities.  The  region  profiles  are 
selected,  within  physical  limits,  to  reproduce  observed  radiances  and/or 
incorporate  other,  related  measurements.  The  utility  of  the  approach  is 
realized  when  systems  studies  address  issues  like  filter  bandpasses  and 
multiple  LOS's.  Signal  and  clutter  intensities  can  be  calculated  for 
different  bandpasses  that  include  the  same  radiators.  Non-linear 
effects  like  those  due  to  opaque  spectral  lines  (e.g.,  clutter  at  60  km 
is  not  relevant  for  limb  viewing  at  4.25  /jm>  are  naturally  included. 
Atmospheric  effects  on  moving  and/or  multiple  targets  can  be  calculated 
as  different  sensor  LOS's  will  traverse  different  parts  of  the  various 
atmospheric  regions. 

The  large-scale  variations  include  diurnal  and  other 
latitude/longitude  effects,  as  well  as  correlations  with  more  global 
parameters  such  as  solar  and  geomagnetic  activity.  Model  atmospheres 
such  as  MSIS,22  and  Jacchia,^^  or  Kenesha^^  empirically  relate 
temperature,  density,  and  major  species  concentrations  to  geophysical 
parameters,  geographic  location  and  universal  time.  Lower  altitudes  and 
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important  IR  species  like  CO2,  O3,  and  NO  would  have  to  be  added.  Their 
densities  could  be  inferred  from  the  major  species  and  temperature 
profiles  using  photochemical  or  empirical  models.  For  example,  O3 
concentrations  above  90  km  may  be  computed  from  the  major  species  via 
known  chemical  reactions. ^>25, 26  msIS  provides  profiles  for  atomic  O, 
which  in  turn  directly  effects  O3  (see  Section  4).  Diurnal  effects  are 
obtained  by  including  photochemical  reactions  in  the  chemical  kinetics. 
The  CO2  profile  depends  more  on  transport  than  chemistry;  assuming  it 
remains  uniformly  mixed,  its  variability  can  be  estimated  by  scaling  to 
the  MSIS  profiles  for  Ar. 

Short-term,  wave-like  responses  of  the  atmosphere  can  be  considered 
as  localized  heating  events  superimposed  on  the  smooth  diffusive 
profiles  of  the  long-term  variations.  These  structures  have  a  wide 
range  of  scale  sizes  and  amplitudes  and  are  characterized  by  irregular 
features  resulting  from  the  dissipative  breakdown  of  waves  into 
turbulent  eddies.  Their  principal  sources  are: 27 

•  turbulence  or  random  fluctuations  in  the  velocity  and  density 
which  serve  to  mix  the  atmosphere, 

•  gravity  wave  structures  that  originate  in  troposphere, 
propagate  through  the  stratosphere  and  ionosphere  into  the 
upper  atmosphere, 

•  diurnal  and  semidiurnal  tidal  motions  that  propagate  from  the 
upper  atmosphere  into  the  middle  atmosphere,  and 

•  winds  and  circulation. 

All  of  these  phenomena  are  closely  related.  For  example,  gravity  waves 
are  associated  with  irregular  wind  components  around  the  mesopause  (70 
to  120  km);  thermospheric  winds  are  the  result  of  thermally  created 
tidal  waves;  and  wind  shears  may  be  a  primary  turbulence  source. 
Geostropic  wind  and  tidal  motion  would  be  considered  as  part  of  the 
average  background  because  their  spatial  and  temporal  scales  are  larger 
than  any  LOS  path. 

Most  atmospheric  flows  are  turbulent  and  occur  on  length  and  time 
scales  that  range  from  several  thousand  kilometers  to  sub  millimeter  and 
years  to  milliseconds,  respectively.  Turbulence  cannot  maintain  itself 
indefinitely,  so  it  must  obtain  energy  from  its  environment.  Turbulence 
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mixes  atmosphere  species  and  opposes  gravity  which  tends  to  separate 
species  vertically  according  to  their  molecular  masses. 

If  not  overwhelmed  by  convection,  buoyant  forces  tend  to  restore 
displaced  parcels  of  air  from  vertical  displacements.  The  air  parcel 
tends  to  oscillate  around  its  undisturbed  position  at  a  frequency  known 
as  the  Brunt -Vaisala  frec[uency.  Typical  periods  are  in  the  vicinity  of 
five  minutes.  Pressure  waves  generated  with  larger  frequencies 
propagate  as  sound  waves,  and  those  with  smaller  frequencies  propagate 
as  internal  gravity  waves.  These  waves  are  excited  near  the  earth's 
surface,  probably  in  wind  systems,  and  the  energy  flow  is  upward, 
reaching  well  into  the  ionosphere.  Their  amplitudes  increase  with 
altitude  and  above  60  km  become  so  large  that  they  dominate  the  wind 
profiles.  Above  about  100  km,  dissipative  mechanisms  like  viscosity, 
thermal  conductivity,  and  eddy  transport  of  heat  and  momentum  break  up 
these  waves  which,  in  turn,  provide  an  energy  source  for  the  upper 
atmosphere.  These  phenomena  have  been  incorporated  into  Global  Climatic 
Models  (GCM)  by  parameterizing  gravity  wave  drag  and  dif fusion. 28“34 

Atmospheric  tides  are  internal  gravity  waves  of  particular 
frequencies.  Tidal  phenomena  usually  imply  gravitational  attraction  by 
the  moon  or  sun,  but  it  can  also  refer  to  waves  caused  by  solar  heating. 
Barometric  pressures  in  the  tropical  latitudes  normally  have  two  maxima 
and  two  minima  per  day  with  a  period  of  12  hours.  At  higher  latitudes 
the  amplitudes  are  smaller  and  are  frequently  dominated  by  noise  from 
passing  weather  systems.  Tidal  wind  patterns  in  the  upper  atmosphere 
also  generate  ionospheric  electrical  currents,  which,  in  turn,  produce 
diurnal  variations  in  the  geomagnetic  field.  These  effects  are  included 
in  some  GCM '  s ,  ,  31 , 35 , 36  g^id  an  initial  approach  for  SHARC  could  use 
different  profiles  with  its  multiple  regions. 

Small-scale  fluctuations  in  the  SHARC  local-regions  model  must  be 
determined  initially  from  field  data.  Fluctuations  in  atmospheric 
density,  composition  and  temperature  provide  the  basis  for  calculating 
the  resultant  IR  clutter.  They  can  either  be  measured  or  obtained  from 
some  first-principles  model  (for  example,  one  for  gravity  waves).  SHARC 
would  calculate  the  excited  state  chemistry  within  each  regions  and  then 
the  LOS  radiances  for  arbitrary  views  and  spectral  bandpasses.  Good 
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agreement  with  observed  radiance  data  would  demonstrate  that  chosen 
power  spectral  density  functions  (PSD)  for  atmospheric  composition  are 
consistent  with  the  data,  but  by  no  means  will  they  form  a  unique  set. 
These  PSD's  then  provide  the  basis  for  predicting  radiances  for  other 
scenarios.  As  an  example  of  this  non-uniqueness,  consider  clutter  in 
the  NO  band  which  might  be  due  to  fluctuations  in  [NO),  [O],  or 
temperature.  If  these  fluctuations  are  in  [O]  or  temperature,  the 
C02(i'2)  band,  for  example,  would  likewise  be  affected  (or  corr-elated )  ; 
if  the  fluctuations  are  mainly  due  to  [NO],  it  would  be  uncorrelated. 
The  crucial  point  is  that  the  chemistry,  radiation  transport,  and 
multiple  region  capability  in  SHARC  provides  a  framework  for  studying 
sources  of  and  models  for  IR  clutter. 
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3. 


SHARC  AURORAL  MODEL 


3 . 1  Model  Deacription 

The  scientific  approach  used  in  SHARC  auroral  model  is 
substantially  based  on  the  Geophysics  Laboratory  Infrared  Spectral 
Radiance  Code  for  the  Auroral  Thermosphere  (AARC).^^  The  calculational 
procedure  for  the  AURORA  module  is  shown  schematically  in  Fig.  8.  The 
primary  function  of  this  module  is  to  perform  a  time  dependent  chemical 
kinetics  calculation  to  determine  auroral  species  number  densities  as  a 
function  of  altitude.  For  C02(i'3)/  the  chemical  kinetics  are  coupled  to 
a  simplified  radiative  transfer  scheme^^  to  obtain  the  excited  state 
number  densities.  The  Gear  method  for  stiff  differential  equations  is 
used  to  perform  the  integration.  This  package  was  supplied  as  part  of 
the  Sandia  CHEMKIN  code,^^  and  thus  it  has  been  thoroughly  tested.  In 
addition  to  the  specification  of  the  auroral  electron  spectra,  the 
AURORA  module  uses  the  input  from  the  "linking"  file  (species 
identification,  stoichiometric  coefficients,  and  rate  constant 
parameters)  to  construct  the  time  dependent  production  rates  for  the 
auroral  species.  After  completing  the  numerical  integration,  the  AURORA 
module  returns  the  vibrationally  excited  number  densities  to  the  SHARC 
main  program.  The  cunbient  and  auroral  populations  are  then  stored  for 


Figure  8.  Major  Components  of  the  Auroral  Model. 
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each  layer  along  the  LOS.  SPCRAD  uses  those  populations  to  calculate 
the  path  radiance  and  transmittance. 

The  remainder  of  this  section  contains  a  general  discussion  of  the 
approach  used  in  the  auroral  model,  a  comparison  of  SHARC  with  the  AARC 
model,  and  an  illustrative  auroral  calculation. 

Calculating  the  radiance  for  a  LOS  through  an  aurora  requires  the 
following: 

•  definition  of  the  strength  of  the  aurora, 

•  specification  of  the  chemical  kinetics  mechanism  for  the 
selected  emitter,  including  primary  and  secondary  electron 
processes,  and 

•  identification  of  the  radiative  transition  and  the 
associated  spectroscopic  data. 

In  addition  to  the  above  information,  a  description  of  the  ambient 
atmosphere  must  also  be  included.  Each  of  the  above  topics  are 
described  in  more  detail  in  the  following  subsections. 

3.1.1  Aurora  Definition 

The  primary  electron  energy  deposition  model  included  in  SHARC  is 
based  upon  work  by  Griin,  Rees,  and  Strickland.  The  model  has  been 

described  in  detail  in  the  AARC  user’s  manual. Routines  describing 
the  energy  deposition  have  been  written  for  SHARC.  These  routines  were 
taken  from  the  AARC  model,  and  modified  to  be  consistent  with  the  SHARC 
structure. 

As  an  example  of  the  primary  electron  energy  deposition  model,  the 
total  ion  pair  production  rate  for  a  class  III"''  aurora  is  shown  in  Fig. 
9.  The  parameters  for  generating  this  ion  pair  production  rate  are 
contained  in  the  AURORA  module  as  part  of  the  data  base.  Of  course,  the 
user  may  choose  to  define  his/her  own  aurora  by  overriding  the  default 
inputs . 
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Figure  9.  Ion  Pair  Production  Rate  for  a  Class  III'*’  Aurora. 


3.1.2  Auroral  Kinetics 

The  auroral  kinetic  model  solves  the  time/energy  dependent  rate 
equations  to  calculate  the  secondary  electron  distributions  and 
reactive/energy  transfer  processes,  and  the  subsequent  enhancement  to 
various  CO2,  NO,  and  N0+  vibrational  states.  The  chemical  reactions  and 
energy  transfer  processes  resulting  from  primary  electron  deposition  in 
the  atmosphere  have  been  described  in  some  detail  in  the  AARC^^  and 
SHARC^  manuals.  The  following  subsections  briefly  describe  the  auroral 
kinetic  mechanisms  employed  by  SHARC,  and  present  comparisons  of  SHARC 
and  AARC  volume  emission  rates  for  the  auroral  emitters. 
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Secondary  Electron  Kinetics 


Since  IR  emiaaions  reaulting  from  direct  electron  impact 
excitation,  as  well  as  recombination  reactions,  depend  on  the  low  energy 
secondary  electron  spectra,  it  is  important  that  the  histogram  treatment 
used  in  SHARC  retain  the  major  features  of  the  secondary  electron 
spectra.  Secondary  electrons  are  formed  via  collisional  ionization  of 
M(=  N2,  O2,  or  O) 

ep  +  M  -*  ep  +  Sfl  products  .  (6) 

From  an  assumed  initial  secondary  electron  energy  distribution,  the 
secondary  electrons  undergo  further  collisions  with  major  species 

eg  +  M  -*  e^  +  M*  (7) 

and 

eg  +  M  -*  +  e4’  +  (8) 

which  result  in  the  electronic  or  vibrational  relaxation  of  the 
secondary  electron  spectra. 

Given  the  ion  pair  production  rate  as  a  function  of  altitude, 
ionization  efficiences,  electronic  and  vibrational  excitation  rate 
constants,  the  secondary  electron  distribution  resulting  from  ionization 
of  the  major  atmospheric  constituents  can  be  calculated.  An  example  of 
a  secondary  distribution  calculated  as  an  input  to  the  AURORA  module  is 
shown  in  Fig.  10.  The  most  interesting  feature  of  this  distribution  is 
the  large  dip  around  3  eV,  which  is  due  to  vibrational  excitation  of  N2 
by  the  secondary  electrons.  The  features  and  altitude  dependence  of 
secondary  electron  distributions  have  been  discussed  in  great  detail 
elsewhere. 

NO  5.3  and  2.7  urn 

Nitric  oxide  in  the  upper  atmosphere  both  during  quiescent  and 
disturbed  conditions  forms  a  major  source  of  infrared  radiation^^  the 
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Figure  10.  Steady  state  Secondary  Electron  Distribution  at  an 
Altitude  of  100  km  for  a  Class  III*  Aurora. 


2.7  to  3.3  fjm  and  the  5.3  to  6  fjm  bands.  When  electrons  interact  with 
the  ambient  atmosphere,  ionization  forms  odd  nitrogen  in  the  form  of 
N(2p),  N(2d),  and  N(^S).  Collisions  between  these  species  and  O2  create 
NO  in  the  upper  atmosphere.  This  production  rate  is  in  turn  found  to  be 
highly  dependent  on  the  ratio  of  N(2d)  to  N(‘^S)  formed  in  the 
interaction  of  electrons  and  N2. 

Nitric  Oxide  is  primarily  produced  by  the  reaction  of  N(2d)  with  O2 
N(2d)  +  O2  -♦  NO(v)  +0  (9) 

and  destroyed  by  reaction  with  ground  state  N(^S) 

N(4s)  +  NO(v)  -»  N2  -*•  O  .  (10) 
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The  dominant  production  mechanisms  for  excited  nitrogen  atoms  are: 
direct  dissociation  of  N2  by  the  primary  electrons  produced  in  the 
electron  beam  (11) >  ion  recombination  reactions  with  secondary  electrons 
(12),  and  charge  transfer  reactions  of  the  primary  ionic  products  (13) 


ep  +  N2  -»  N(4s)  +  N(2d)  -t-  ep  (11a) 

ep  +  N2  -•  N(2d)  +  N'^  +  eg  +  ep  (Hb) 

NO"^  +  eg  N(2d)  +  O  (12) 

N2'*‘  +  0-»N(2d)+N0.  (13) 


Nitrogen  atoms  are  also  subject  to  quenching  by  other  species  which 
compete  with  the  reactions  of  Eqs.  (9)  and  (10)  in  the  production  and 
destruction  of  NO 


N(2d)  +  N2 

-  N(4s)  +  N2 

(14) 

N(2d)  +  0 

-*  N(4s)  +  0 

(15) 

The  vibrational ly  hot  N0(v)  can  be  quenched  by  O2 

NO(v)  +  O2  -»  NO(v-l)  +  O2  .  (16) 

Similarly,  the  upper  levels  can  emit  radiation,  and  thus  cascade  down 
the  vibrational  ladder. 


NO(V) 

-*  NO(v-l)  +  hi/ 

(17a) 

NO(v) 

-*  NO(v-2)  +  hi/ 

(17b) 

Radiative  decay  rates  from  the  high  v  states  by  Av=2  are  comparable  to 
that  for  Av=l. 


Laboratory  and  field  meaeurement  programs  have  indicated  the 
existence  of  extremely  high  rotational  excitation  (rotational 
temperatures  on  the  order  of  10000  K)  for  vibrational  states  up  to  v=5. 
There  is  some  evidence  that  the  mechanism  for  the  formation  of  this 
rotationally  hot  NO  is  from  the  reaction 

N(2p)  +  O2  -*  N0t{v)  +  O  (18) 
although  this  process  is  not  currently  treated  in  SHARC. 

CO2  4.3  urn 

There  are  two  ways  to  excite  CO2  efficiently  in  an  aurorally 
excited  atmosphere:  direct  t  .ctron  excitation,  and  the  transfer  of 
energy  from  vibrationally  excited  Nj  to  CO2 .  Different  models  of  N2 
excitation  due  to  secondaries  have  as  much  as  a  factor  of  three 
disagreement.  Assumptions  differ  on  the  precise  shape  of  the  secondary 
electron  energy  spectrum  and  mechanisms  assumed  for  thermal izing  these 
electrons. 

In  order  to  calculate  the  intensities  of  the  prompt  and  delayed 
emissions,  several  assumptions  are  made.  Important  parameters  are  the 
ion  pair  production  rate  and  the  aurora  duration.  The  total  amount  of 
energy  put  into  the  atmosphere  determines  the  amount  of  energy  stored  in 
N2,  while  the  ion  pair  production  rate  determines  the  number  of 
secondary  electrons  present  to  transfer  energy  directly  to  CO2 .  These 
reactions  are  represented  in  SHARC  by: 


Ss  +  N2 

-*  eg  +  N2(v) 

(19a) 

N(4s)  + 

NO 

-»  0  +  N2  ( v) 

(19b) 

0(1d)  + 

N2 

-*  0  +  N2  ( V) 

(19c) 

N2(v)  + 

CO2 

-*  N 2  C02(t'3) 

(20) 
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By  binning  the  secondary  electrons  eg  and  e^,  both  the  energy  spectra 
and  the  N2(v)  are  simultaneously  treated.  The  number  of  vibrationally 
hot  N2  molecules  created  in  this  way  was  found  to  vary  with  altitude  and 


dose. 


The  mechanisms  for  the  prompt  emission  at  4.3  pm  are  given  by 
®s  +  CO2  -*  e^  +  C02(V3)  (21) 

C02(i'3)  -*  CO2  +  hi/  .  (22) 

NO'*'  4.3  urn 

NO'*'  emission  in  the  4.3  pm  region  has  not  been  uneunbiguously 
identified  in  the  region  to  the  red  of  CO2  1/3  due  to  the  lack  of 
high-sensitivity,  high  resolution  spectral  measurements.  Previous 
models  have  assumed  the  primary  source  of  NO'*  formation  is  from  the  ion 
molecule  reaction  of  N'*  with  03.  The  primary  reaction  set  used  in  these 
calculations  is  as  follows: 


ep  +  N2  N-*  +  jN(2d)  +  jN(4s)  +  eg  +  ej 


N'*'  +  O2  -»  NO'*'(v)  +  O 


The  formation  of  NO'*  competes  with  the  charge  exchange  reaction. 


N"*  +  O2  -»  Oj  +  N 


The  vibrationally  excited  NO'''(v)  can  be  quenched  by  N2 


NO'*'(v)  +  N2  -*  NO'*'(v-l)  <■  N2 


as  well  as  emit  radiation,  and  thus  cascade  down  the  vibrational  ladder 


NO'*'(v)  -»  N0'''(v-1)  +  hi/ 
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3.1.3  Radiative  Transfer  Model  for  C02(4.3  ^Jm) 

The  emission  from  NO  and  NO"*"  is  both  optically  thin  and  prompt. 
Thus,  coupling  of  the  radiative  transfer  and  collisional  processes,  is 
not  required,  and  the  solutions  to  the  mechanisms  given  previously 
suffice  in  the  calculations.  However,  for  002(4.3  pm),  emission  below 
100  km  is  both  optically  thick  and  delayed  due  to  the  N2-C02(i^3) 
vibrational  energy  transfer.  Thus,  to  treat  the  radiative  transfer 
correctly  requires  the  coupling  of  the  CO2  time-dependent  chemical 
kinetics  with  the  radiative  transfer  scheme.  Unfortunately,  such  a 
treatment  in  SHARC  is  prohibitively  time  consuming  and  is  not  possible 
at  this  time. 

Kumer^^  has  presented  a  model  in  which  he  has  shown  how  to  uncouple 
the  radiative  transport  from  the  time-dependent  kinetic  processes,  and 
which  he  has  termed  the  approximation.  The  approximation  has  also 
been  discussed  in  some  detail  in  the  AARC  report.  In  this 

approximation,  the  chemical  kinetic  equation  for  CO2  becomes 

^(C02(001)lj.  -  ke(e3]i[C02{000))i  +  ki  (Nj  ]  ilCOj  (000)  )  i 

-  ki(N2li(CO2(001) -  kM(M)i(CO2(001) )i  (28) 

-  (l-2wji)A4  3(002(001)  )i  +  ka[CO2(000)  )i 

j 

where  the  first  four  terms  in  Eg.  (28)  result  from  the  direct  secondary 
electron  excitation  of  CO2 #  the  vibrational  energy  transfer  of  CO2  with 
N2 ,  and  the  quenching  of  002(001)  by  the  atmosphere.  The  last  two  terms 
represent  the  effects  of  radiation  transport  on  the  chemical  kinetics. 
The  first  effect  (given  by  the  fifth  term)  is  simply  a  modification  of 
the  Einstein  coefficient,  3,  due  to  the  fact  that  not  all  photons 
emitted  from  the  i^^  layer  escape  the  atmosphere.  The  second  effect  is 
included  as  an  external  source,  represented  by  the  rate  constant  k^,  and 
accounts  for  the  excitation  of  CO2  in  the  i*^^  layer  by  radiation  emitted 
from  all  other  layers. 
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An  analogous  rate  equation  for  the  excitation  of  N2  is 


=  ke[esli(N2]i  +  ( N2  1  i  1^02  (001 ))  i 

-  >tl(N2]i(CO2(000)  li  +  k2(N2)itCO2(011)  li 

-  Jt2lN2)i[C02(010)]i  -  kMlMJilNjJi  (29) 

+  k3e(N(4s) )i[N0]i  +  k4e(0(lD) ]i(N2]i 

where  the  first  six  terms  in  Eq.  (29)  result  from  the  direct  secondary 
electron  excitation  of  N2/  the  vibrational  energy  transfer  of  CO2(001) 
and  CO2(011)  with  N2,  and  the  quenching  of  by  the  atmosphere.  The 
last  two  terms  in  Eq.  (29)  represent  the  effective  rate  of  formation  of 
N2  from  aurorally  produced  species.  The  relative  importance  of  these 
terms  in  producing  is  discussed  in  the  next  section. 

3 . 2  Coinpariso..  SHARC  with  AARC 

This  section  presents  a  brief  validation  of  the  chemical  kinetics 
mechanisms  used  by  SHARC  by  comparing  the  efficiency  of  producing 

photons  per  ion  pair  predicted  by  SHARC  with  those  from  AARC.l^  The 

intent  of  this  section  is  not  to  discuss  the  particular  features  of  the 
efficiencies  since  this  has  been  done  in  some  detail  in  the  AARC  report. 
Efficiencies  are  computed  by  SHARC  for  NO(Av=l),  N2(v=l),  and  NO'''(Av=l) 
assuming  an  IBC  III  aurora  and  a  standard  atmosphere  which  is  included 
with  SHARC. 

The  altitude  and  time  dependence  of  the  NO(Av=l)  photons/ion  pair 
predicted  by  SHARC  is  compared  with  the  steady  state  AARC  prediction  in 
Fig.  11.  At  altitudes  less  than  115  km,  SHARC  indicates  that  the 
efficiency  is  indeed  time  independent  and  in  excellent  agreement  with 
the  AARC  calculation.  However,  at  altitudes  greater  than  115  km,  the 
number  of  photons  produced  per  ion  pair  shows  a  significant  time 
dependence  and  is  almost  a  factor  of  2  greater  than  the  AARC  efficiency. 
This  effect  results  almost  entirely  from  the  collisional  excitation  of 
ground  state  NO  (which  is  continually  produced  by  the  aurora)  by  O 
atoms.  However,  in  AARC,  an  initial,  time  independent  NO  number  density 
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profile  is  assumed,  and  the  steady  state  NO  vibrational  populations  are 
then  normalized  to  the  input  profile.  Although  there  is  a  sizeable 
difference  between  SHARC  and  AARC  at  high  altitudes,  the  effect  of  this 
difference  is  relatively  unimportant  for  most  LOS  radiances  because 
emissions  from  the  ambient  atmosphere  are  dominant  at  these  altitudes. 

One  significant  difference  between  SHARC  and  AARC  is  the  rate 
constant  assumed  for  the  process 

ep  +  N2  -*  N(4s)  +  N(2d)  -t-  ep  (11) 

which  results  in  almost  a  factor  of  two  difference  in  the  predicted 
efficiencies  as  shown  in  Fig.  12.  The  SHARC  value  for  Eg.  (11)  is  based 
upon  work  done  in  the  GL  EXCEDE  program  and  is  more  than  a  factor  of  two 
less  than  the  value  assumed  by  AARC. 


EPP i c i ency 


(photons/ i on  -po i r) 


Figure  11.  Comparison  of  the  SHARC  NO(Av=l)  Time-Dependent 
Photons/Ion  Pair  with  the  AARC  Efficiency. 
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Figure  12.  Comparison  of  the  SHARC  NO(Av=l)  Predictions  of 
Photon/Ion  Pair  Efficiencies  Using  the  EXCEDE  and  AARC 
Rate  Constants  for  Eg.  (11). 


The  dominant  mechanism  for  the  CO2  4.3  pm  emission  is  the  transfer 
of  energy  from  vibrationally  excited  N2  to  CO2.  The  most  important 
processes  for  formation  of  are  the  direct  secondary  electron 
excitation  of  N2»  the  reaction  of  N  atoms  with  NO,  and  electronic  to 
vibrational  energy  transfer  between  0(1d)  and  N2  (see  Eqs.  (19)).  The 
altitude  dependence  of  the  SHARC  efficiency  due  to  each  of  these 
processes  is  shown  in  Fig.  13.  For  the  assumed  mechanisms  and  rate 
constants,  it  can  be  seen  that  the  direct  electron  excitation  of  N2  is 
of  comparable  importance  to  the  chemical  production  mechanisms. 
Furthermore,  there  is  a  only  a  slight  altitude  dependence  to  the 
efficiencies.  The  AARC  prediction  of  the  N2  efficiencies  for  the  same 
processes  are  shown  in  Fig.  14.  For  the  chemical  mechanisms  (19b  and 
19c),  the  agreement  between  SHARC  and  AARC  is  no  worse  than  30%, 
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Figure  13.  SHARC  Prediction  of  the  N2(v)  Efficiencies  (Quanta/Ion 
Pair)  for  Several  Production  Mechanisms. 
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Figure  14.  AARC  Prediction  of  the  N2(v)  Efficiencies  (Quanta/Ion 
Pair)  for  Several  Production  Mechanisms. 


EFPiciency  (photons/ i on -po i r) 

Figure  15.  Comparison  of  the  SHARC  and  AARC  NO''’(Av=l) 
Efficiencies  (Photons/Ion  Pair). 


although  the  detailed  profile  of  the  efficiency  due  to  N(^S)+NO  is  quite 
different  due  to  the  manner  in  which  the  NO  profile  is  generated  during 
an  aurora.  The  most  substantial  difference  between  SHARC  and  AARC  is 
the  secondary  electron  excitation  efficiency  at  altitudes  greater  than 
110  km,  which  is  not  understood  at  the  present  time. 

The  final  comparison  of  SHARC  and  AARC  efficiencies  is  shown  in 
Fig.  15  for  NO'*'(Av=l).  The  agreement  between  SHARC  and  AARC  is 
excellent  over  the  entire  altitude  regime. 


3 . 3  Illustrative  Calculation 


As  an  example  of  the  auroral  enhancement  of  NO  5.3  fjtn  radiation,  a 
calculation  has  been  performed  for  an  IBC  II  aurora,  appropriate  to  the 
FWi19  experiment.  The  calculation  uses  an  augmented  HITRAN  data  file,!'^ 
which  has  been  expanded  to  include  lines  for  high  NO  vibrational  states 
and  NO*. 
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Volumetric  Emission  Rote  (photons/s/cm^/sr) 


Figure  16.  SHARC  Prediction  of  Ambient  and  Auroral  IBC  II 
Altitude  Dependent  Volumetric  Emission  Rates  for 
NO(Av*l ) . 


The  relative  contributions  of  the  ambient  and  auroral  components  to 
the  volumetric  emission  rate  as  a  function  of  altitude  are  illustrated 
in  Fig.  16  for  NO(Av=l).  The  ambient  emission  increases  with  altitude 
until  approximately  125  km  and  then  begins  to  decrease.  This  behavior 
is  due  to  O+NO(v=0)  collisional  excitation,  which  depends  on  the  O  and 
NO  number  density  profiles  and  the  ambient  temperature  (see  Eq.  5).  The 
altitude  dependence  of  the  auroral  emission  is  dictated  by  the  N(^D) 
altitude  profile,  which  in  turn  is  determined  by  the  ion  pair  production 
rate  ( see  Fig.  9) . 

A  comparison  of  a  SHARC  auroral  calculation  to  FWI  data^^  is  shown 
in  Fig.  17.  The  data  were  taken  from  a  sounding  rocket  launched  from 
Poker  Flat,  Alaska  in  April  1973  under  auroral  conditions.  These 
measurements  provided  the  first  high  spectral  spectral  resolution  IR 
data  for  an  aurorally  disturbed  atmosphere.  Ground  measurements  prior 
to  launch  indicated  an  auroral  strength  of  50  kiloRayleighs,  but  it  had 
weakened  to  less  than  12  kR  by  launch  time.  The  spectral  data  have  a 
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resolution  of  about  1  cm“l.  The  model  calculations  are  for  a  Class  II 
aurora  of  10  kR.  Although  the  calculations  fall  about  a  factor  of  four 
below  the  data,  the  results  are  quite  promising  for  an  ab  initio 
calculation  that  starts  with  electron  dosing  of  the  ambient  atmosphere. 
During  auroral  events,  the  NO  and  N0+  populations  increase 
significantly,  and  we  ascribe  the  difference  to  strong  predosing  and 
population  enhancement  prior  to  the  measurements,  which  are  not  modeled 
in  the  present  version  of  SHARC.  It  is  also  possible  that  the 
uncertainty  in  the  rate  of  N(2d)  creation  from  Eg.  (11a)  as  shown  in 
Fig.  12  is  responsible  for  part  of  the  discrepancy.  Further  detailed 
analysis  of  the  full  data  set  is  required  to  clarify  this  point. 
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Figure  17. 


Comparison  of  the  SHARC  Spectral  Radiance  with  FWI 
Data  for  a  Near  Zenith  Line-of -Sight . 
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4.  ANALYSIS  OF  SPIRIT  1  DATA 


The  Spectral  Infrared  Interferometric  Telescope  (SPIRIT  1) 
experiment,  launched  in  1986  from  Poker  Flat,  Alaska,  generated  a  large 
and  valuable  body  of  earth  limb  LWIR  spectral  data  during  auroral 
conditions . Emissions  were  observed  from  CO2  ^2'  fundamental, 
and  the  O3  t/3  fundamental  and  hot  bands.  These  data  vBre  analyzed  by 
SSI  as  part  of  an  overall  effort  to  improve  our  understanding  of  LWIR 
backgrounds  and  to  validate  model  calculations. 

Data  analysis  under  this  effort  included  data  compilation  and 
selection,  identification  of  background  contamination  sources, 
contamination  subtraction,  spectral  modeling  using  atmospheric  radiation 
codes,  limb  inversions  using  a  specially-developed  least-squares  method, 
temperature  profile  determinations  using  a  novel  moment  method,  and 
density  profile  determinations  using  photochemical  models.  In  the 
course  of  this  work,  several  upgrades  to  the  raw  data  reduction  scheme 
developed  by  Boston  College  were  also  implemented.  These  included  the 
determination  and  application  of  wavelength  response  curves  and  scale 
factors  for  absolute  radiance  calibration,  and  tangent  height 
corrections  based  on  data  from  the  on-board  horizon  and  star  sensors. 
The  results  of  the  SPIRIT  1  data  analysis  program  are  summarized  in  four 
publications,  which  are  described  briefly  in  the  following  subsections. 

A  general  conclusion  to  be  drawn  from  the  SPIRIT  1  observations  is 
that  the  LWIR  limb  radiances  as  well  as  the  average  densities  of  the 
pertinent  species  (O,  O3,  CO2,  and  NO)  in  this  experiment  were  typically 
within  a  factor  of  two  of  current  non-auroral  model  predictions. 
Features  such  as  NO  hot  bands  that  are  known  to  arise  from  auroral 
excitation  channels  were  not  observed  in  SPIRIT  1  due  to  their  weakness 
relative  to  the  quiescent  atmospheric  limb  background,  or  due  to 
instrumental  effects  such  as  the  presence  of  strong  CO2  ^2  radiation  or 
background  contamination.  However,  evidence  may  exist  for  more  subtle 
effects,  such  as  O  and  H  atom  depletion,  associated  with  atmospheric 
disturbances.  It  is  suggested  that  future  auroral  observation 
experiments  utilize  nadir-  or  slant-viewing  rather  than  limb-viewing 
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geometries  to  minimize  off-axis  contamination  as  well  as  the  normal 
atmospheric  background  radiation. 

An  important  result  from  SPIRIT  1  is  the  acquisition  of  high- 
quality,  well-resolved  ozone  hot  band  spectra  in  the  10-13  pm 
atmospheric  window.  Analysis  of  these  spectra  have  resulted  in  new  hot 
band  center  assignments,  as  well  as  a  quantification  of  the  relationship 
between  the  hot  band  intensity  and  the  O-atom  recombination  rate. 

4.1  9-12  um  Ozone  Emission 

Spectra  of  ozone  1/3  emission  observed  in  the  auroral,  nighttime 
atmosphere  during  the  SPIRIT  1  experiment  were  analyzed.^®  The  well 
resolved  hot-band  spectra  reveal  spectral  structure  not  previously 
observable,  and  provide  a  critical  test  of  ozone  radiance  models.  The 
limb  radiance  profile  is  in  good  agreement  with  previous  data  below  -95 
km,  but  above  this  altitude  it  falls  off  more  rapidly.  Results  of  this 
study  include  limb  spectra,  limb  radiances  and  volumetric  emissions  in 
various  band  passes.  Also  a  nighttime  ozone  profile  was  derived. 

Examples  of  the  ozone  spectra  are  shown  in  Fig.  18.  The 
(OOl)-(OOO)  (cold)  band  appears  as  the  strong  double-peaked  feature. 
The  1/3  hot  bands,  which  arise  from  O  +  O2  +  M  recombination,  appear  as 
the  tail  to  the  red  of  the  cold  band. 

Figure  19  compares  ozone  density  profiles  obtained  from  several 
different  sources;  a  calculation  based  on  the  SPIRIT  1  cold  band 
radiance  profile  with  an  estimated  Ty;  the  profile  assumed  for  the  HAIRM 
code;H  and  a  photochemical  model  calculation  with  assumed  O  and  H 
profiles.  The  f actor-of-two  consistency  between  these  profiles  up  to 
altitudes  of  100  km  supports  the  claim  that  the  essential  physics  and 
chemistry  of  ozone  formation,  destruction,  and  radiation  at  these 
altitudes  are  reasonably  well  understood.  An  improved  analysis  resulted 
in  more  accurate  O,  H  and  Ty  profiles.'*^ 

The  ozone  hot  band  spectrum  is  briefly  discussed.  The  need  to 
update  the  current  hot  band  center  frequencies  in  HAIRM  and  SHARC  was 
noted,  and  a  modified  formula  for  ozone  vibrational  energy  levels  was 
proposed. 
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18.  Ozone  Limb  Spectra  from  SPIRIT  1  with  Kaiser-Bessel 
Apodization. 
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Figure  19 


Comparison  of  Standard  Nighttime  Ozone  Density 
Profiles  with  Profile  Derived  from  SPIRIT  1  Data. 


4 . 2  Identification  of  O3  High  Vibrational  States 

O3  hot  band  chemiluminescence  spectra  obtained  in  SPIRIT  1  were 
analyzed  in  detail  with  the  aid  of  laboratory  data  and  vibrational 
energy  and  band  contour  predictions.^^  Band  centers  and  energies  for 
previously  unreported  vibrational  states  of  ozone  having  up  to  seven 
quanta  were  derived.  The  assignments  are  supported  by  the  predictions 
of  an  improved  Darling-Dennison-type  energy  level  formula  which  behaves 
realistically  towards  dissociation.  Bands  associated  with  bending  as 
well  as  stretching  excitation  appear  due  to  an  accidental  vibrational 
resonance  and  a  non-mode-specific  population  distribution. 

The  band  centers  and  vibrational  energies  derived  from  this 
analysis  are  summarized  in  Table  2.  Eight  bands  whose  isolated  Q 
branches  were  observed  in  the  SPIRIT  1  spectrum  appear  in  the  "Current" 
column;  their  centers  are  accurate  to  within  -2  cm~^  or  better. 

A  synthetic  spectrum  was  generated  by  least-squares  fitting  to  the 
observed  spectrum  using  eight  independent  population  components.  The 
overall  agreement  between  the  synthetic  and  observed  spectra  (Pig.  20) 
is  quite  good,  considering  the  approximations  in  the  spectral  model. 

4 . 3  C02(i'7)  and  N0fAv=ll  Emissions 

SPIRIT  1  spectral  data  for  CO2  >^2  emissions  were 

analyzed. These  data  cover  auroral  intensities  from  several  kR  to 
over  100  kR  at  391.4  nm,  and  tangent  heights  of  -70-200  km.  Volumetric 
emission  and  kinetic  temperature  profiles  have  been  obtained  using 
linear  least-squares  inversion  procedures.  The  CO2  t^2  radiance  profile 
is  found  to  be  similar  to  previous  observations.  The  NO  5.3  fjm  emission 
is  weaker  than  usual,  which  is  ascribed  mainly  to  a  low  thermospheric 
temperature.  Auroral  enhancement  of  NO  hot  bands  was  not  observable 
with  the  available  sensitivity. 
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Table  2. 

Energies 

&  Band 

Centers 

for  O3  Vibrational  States 

!  STATE 

ENERGY  (cm- 

1) 

BAND 

CENTER 

Calc. 

Obs. 

Current 

Calc. 

Obs. 

Current 

;  003 

3045 

304611® 

3046 

987 

98713*3 

988 

004 

3999 

4000 

954 

95213*3 

954® 

005 

4910*^ 

91ld 

4918® 

491811^ 

4919 

919® 

91613*3 

919 

311 

4902«i 

4894® 

489511^ 

4894 

8959'® 

8949 

006 

5767 

5761110^ 

5771 

841*^ 

849® 

852 

8739-® 

8779 

102 

3085 

308411® 

3083 

975 

97113*3 

972 

103 

4025 

4019llf 

4020 

940 

93913*3 

937® 

104 

4925 

4926 

900 

90413*3 

906® 

105 

5774 

5778 

849 

852® 

013 

3695 

369711® 

3697 

970 

971 

014 

4631 

4633 

936 

936 

015 

5519 

5517 

888 

884® 

016 

6342 

6345 

823 

828® 

112 

3743 

3743 

958 

958 

113 

4663 

465911^ 

4660 

920 

917 

114 

5537 

5536 

874 

876® 

115 

6351 

6353 

814 

817® 

®  Barbe,  et 

al.  (Ref. 

46);  absorption 

^  Rawlins  and  Armstrong  (Ref.  47);  emission 
^  Q  branch  identified  in  SPIRIT  1  spectrum 
^  Excludes  (311)-(005)  resonance 
®  Includes  (311)-(005)  resonance 

^  Contour  analysis  of  Damon  et  al.  (Ref.  48)  absorption  spectrum 
9  Crossover  band 

^  Imre  (Ref.  49);  Raman  spectrum 


Examples  of  the  data  appear  in  Figs.  21  thru  23.  Figure  21  shows 
the  CO2  1^2  radiance  versus  altitude,  and  Fig.  22  shows  representative 
CO2  spectra  at  different  altitudes.  The  NO  radiances  appear  in  Fig.  23. 
Figure  24  shows  the  results  of  a  determination  of  the  mesospheric 
temperature  profile  '  'om  the  NO  band  shape  using  a  new  method  based  on 


40 


TANGENT  HEIGHT  (kn)  ®  REL .  INTENSITY 
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20.  Comparison  of  Observed  and  Synthetic  Unapodized  Ozone 
Spectra. 
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Figure  21.  CO2  *^2  Limb  Radiance  Profile  From  SPIRIT  1,  and 
Average  of  SPIRE  Data. 
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Figure  22.  Typical  CO2  ^2  Spectra  at  Various  Tangent  Heights  with 
Kaiaer-Bessel  Apodization  (11  cm~l  Resolution). 
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Figure  23.  NO  Av=l  (5.3  ^m)  Limb  Radiance  Profile,  Showing 

Second-Degree  Fit  and  Full  Limb-Corrected  SPIRIT  1 
Data.  SPIRE  Full  Limb  Data  are  also  Plotted. 
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Figure  24,  Limb  and  Local  Kinetic  Temperature  Profiles  from 
SPIRIT  1  NO  Av=l  Data.  Error  Bars  Show  Approximate 
Uncertainty  in  Local  Temperature. 

spectral  moments.  To  within  the  rather  large  error  bars  which  arise 
from  scatter  in  the  band  shape,  the  determined  profile  is  in  agreement 
with  the  MSIS  model. 

4.4  O  and  H  Profiles  Derived  for  SPIRIT  1  and  SPIRE 

A  method  was  developed  for  deriving  atmospheric  O  and  H  atom 
density  profiles  in  the  mesopause  region  using  ozone  1^3  emissions. 44  xhe 
10-12  fim  hot  band  intensity  determines  the  O  +  O2  +  M  recombination 
rate,  while  the  nighttime  cold  band  intensity  determines  the  chemical 
destruction  rate  and  hence  the  H  atom  density.  Cold  band  emissions 
under  sunlit  conditions  (in  the  SPIRE  experiment)  have  been  used  to 
provide  an  absolute  calibration.  Derived  density  profiles  are  presented 
for  the  SPIRIT  1  and  SPIRE  rocket  experiments,  and  compare  reasonably 
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with  KSIS-86  model  predictions.  The  somewhat  lower  atom  densities  found 
in  the  auroral  SPIRIT  1  conditions  may  reflect  local  transport  effects. 

The  computed  O  and  H  atom  profiles  for  the  SPIRIT  1  experiment  are 
shown  in  Figs.  2S  and  26.  The  SPIRE  O  atom  density  is  2-3  times  larger, 
in  better  agreement  with  MSIS.  The  estimated  absolute  accuracy  of  the  o 
atom  profile  is  within  a  factor  of  two,  while  the  H  atom  profile  is 
believed  to  be  more  accurate  than  this. 
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Figure  25.  O  Atom  Profiles  Calculated  from  the  Inverted  SPIRIT  1 
Hot  Band  Data. 
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Figure  26.  H  Atom  Profiles  Calculated  from  the  SPIRE  and  SPIRIT  1 
Data.  MSIS-86  Model  Profiles  are  Essentially 
Identical  for  Both  Experiments. 
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5.  ANALYSIS  OF  EXCEDE  DATA 


The  goal  of  EXCEDE  III  is  to  provide  a  data  base  to  reduce 
uncertainties  in  modeling  electron-dosed  atmospheres  this  end  the 
response  of  the  natural,  quiescent  atmosphere  to  dosing  by  energetic 
electrons  is  studied.  A  beam  of  electrons  having  a  nominal  initial 
energy  of  2.6  keV  was  injected  into  the  upper  atmosphere  between  80  and 
115  km.  Optical  instruments  observed  the  interaction  of  these  electrons 
with  the  atmosphere  from  both  on-board  the  electron  accelerator  bearing 
module  and  from  off-board  platforms.  EXCEDE  III  is  the  fifth  flight  in 
the  EXCEDE  program  which  has  demonstrated  the  principle  of  artificially 
creating  a  electron-excited  cloud  in  the  upper  atmosphere. 

Four  Circular  Variable  Filter  Spectrometers  (CVFs)  were  flown  on 
EXCEDE  III;  one  was  located  on  the  accelerator  module  to  look  at 
emissiono  near  the  electron  beam  source,  while  three  were  placed  on  the 
sensor  module.  Of  these  three,  one  CVF  was  pointed  at  the  electron  beam 
deposition  region,  looking  at  the  prompt  emission,  a  second  looked  in 
the  region  8°  behind  the  electron  beam,  corresponding  to  emissions 
nominally  0.2  seconds  in  the  afterglow,  and  a  third  looked  20°  behind 
the  electron  be£un,  measuring  emissions  0.5  seconds  in  the  afterglow. 
The  CVFs  had  continuous  spectral  coverage  from  2.4  to  21  fim  and  3%  FWHM 
spectral  resolution.  Data  from  the  three  sensor  CVFs  measured  emissions 
which  provide  important  information  on  the  time-dependent  chemistry 
responsible  for  IR  radiation.  With  the  aid  of  interferometer 
measurements,  emissions  have  been  positively  identified.  Infrared 
spectra  of  CO2,  CO,  NO,  and  O3  were  seen  both  during  beam-on  and 
beam-off  periods.  Time  histories  of  infrared-active  bands  are  studied. 

Current  uncertainties  in  the  phenomenological  models  based  on 
reaction  rate  coefficients  and  calculated  electron  concentrations  that 
affect  infrared  emissions  are  thought  to  be  a  factor  of  2  or  3  in  the 
best  known  situations  and  much  worse  in  many  other  cases.  Thus,  it  is 
important  to  validate  both  the  emission  rate  and  the  mechanisms  of 
formation  of  important  IR  radiators  in  the  disturbed  atmosphere. 
Particular  species  of  interest  are  C02(^'2  *^3)'  N0(Av  =  1  and  Av  = 
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2),  CO(Av=l),  NO*,  O3,  and  other  oxides  of  nitrogen  such  as  N2O  and  NO2 . 
Unknowns,  or  disputed  reaction  rates,  include  the  quenching  of  precursor 
species  such  as  N(2d)  by  0(2p),  fractional  yields  of  species  such  as 
N(2p)  and  N2(A),  effects  of  these  species  in  the  disturbed  atmosphere, 
and  creation  mechanisms  for  species  such  as  O3  and  N2O.  Emissions  due 
to  very  high  rotation  levels  which  are  not  quenched  in  the  field 
experiment  may  impact  the  atmospheric  bandpasses  being  considered  for 
possible  surveillance  instruments. 

This  section  discusses  a  first-look  analysis  of  the  infrared  data 
obtained  from  the  Circular  Variable  Filter  Spectrometers  (CVFs)  flown  on 
the  EXCEDE  III  sensor  module.  Only  a  small  fraction  of  the  available 
data  and  model  calculations  is  presented  here.  A  post-flight 
calibration  of  the  CVF  instruments  has  recently  been  completed,  and, 
therefore,  the  data  shown  here  is  not  yet  final.  Although  much  of  this 
work  is  preliminary,  qualitative  observations  from  the  data  and  kinetic 
model  will  no  doubt  survive  a  more  rigorous  quantitative  analysis. 

5. 1  Sensor  Module  CVF  Instrumentation 

The  purpose  of  the  CVF  instruments  is  to  obtain  spectral 
information  on  the  prompt  and  delayed  (afterglow)  emission  resulting 
from  the  irradiation  of  the  atmosphere  by  an  electron  beam.  The 
approximate  scan  time  for  each  CVF  was  1.8  sec.  The  spectral  coverage 
was  2.4-21  pm  with  a  3%  resolution. 

As  shown  in  Fig.  27,  the  arrangement  of  three  CVFs  on  the  sensor 
module  at  0°,  8°,  and  20°  with  respect  to  the  electron  beam  allows 
spectral  data  to  be  obtained  not  only  for  prompt  emission,  but  also  for 
the  afterglow  region.  In  particular,  for  the  geometry  of  the  EXCEDE  III 
flight,  spectral  data  is  obtained  for  up  to  0.5  sec  after  the  electron 
dosing. 
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Figure  27.  Sensor  Module  CVFs  Fields-of-View. 

5 . 2  Representative  CVF  Data 

During  the  course  of  the  experiment,  each  of  the  thrae  sensor 
module  CVFs  obtained  approximately  88  spectral  scans.  Two  complete 
spectral  scans  were  made  during  each  beam-on  period  and  one  scan  during 
beam-off.  The  altitude  range  over  which  the  scans  were  taken  was 
approximately  85  km  to  115  km  (apogee).  Although  excellent  quality 
spectra  were  obtained  for  both  electron-beam  initiated  emission  and 
emission  from  the  quiescent  atmosphere,  it  should  be  realized  that  these 
data  are  still  preliminary.  In  particular,  no  data  are  shown  from  the 
0°  CVF  due  to  a  calibration  problem.  The  CVFs  have  only  recently 
undergone  a  post-flight  calibration,  and  the  recalibrated  data  are  not 
yet  available. 

The  effect  of  the  electron-beam  dosing  in  the  2. 4-6.0  pm  is  shown 
in  Fig.  28.  The  data  are  from  the  8°  CVF  looking  approximately  0.2  sec 
into  the  afterglow.  The  time  period  shown  includes  three  electron  gun 
pulses  corresponding  to  an  altitude  range  of  113  to  107  )tm.  The 
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Figure  28.  Electron  Beam  Enhancement  of  NO(Av=l)  and  N0(Av=2) 
Radiance  for  Three  Electron  Gun  Pulses. 

enhancement  of  the  NO(Av=l)  and  NO(Av=2)  emission  at  5.3  /jm  and  2.7  pm 
is  clearly  shown.  Although  not  evident  from  the  figure,  there  is 
significant  hot  band  emission  at  wavelengths  greater  than  5.3  pm  (from 
vibrational  states  >  2)  which  is  not  present  in  the  ambient  emission. 

Figure  29  illustrates  the  effect  of  the  electron-beam  dosing  in  the 
8.5-10.5  pm  spectral  region.  Again,  the  data  are  from  the  8°  CVF.  This 
time  period  includes  four  electron  gun  pulses  corresponding  to  an 
altitude  range  of  110  to  99  km.  As  expected,  the  O3  (9.6  pm)  intensity 
becomes  significant  below  approximately  105  km.  In  contrast  to  NO, 
there  does  not  appear  to  be  any  significant  enhancement  of  the  O3 
(9.6  pm)  emission  over  the  ambient  background.  However,  this  does  not 
mean  there  is  no  O3  (9.6  pm)  enhancement  due  ’_o  the  electron  beam,  only 
that  the  background  is  large  enough  to  mask  any  electron  beam 
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Figure  29.  Effect  of  Electron  Beam  on  O3  (9.6  ^Jm)  Radiance  During 
Four  Electron  Gun  Pulses. 

excitation.  In  fact,  recent  laboratory  experiments  at  the  Geophysics 
Laboratory  have  indicated  that  O3  (9.6  pm)  is  very  efficiently  excited 
by  electron  bombardment.^*^ 

The  CVF  data  also  provide  a  wealth  of  information  on  the  time  and 
altitude  dependence  of  quenching  and  vibrational-to-vibrational  energy 
transfer.  For  example,  Fig.  30  shows  the  effects  of  the  time-dependent 
collisional  processes  by  comparing  the  8°  and  20°  CVF  data.  Recall  that 
the  8°  and  20°  CVFs  look  approximately  0.2  and  0.5  sec,  respectively, 
after  beam  excitation.  The  effect  of  quenching  on  the  prompt  emitters, 
such  as  NO(Av=l  and  2)  at  5.3  and  2.7  pm  and  CO(Av=l)  at  4.7  pm  is  quite 
dramatic.  However,  the  behavior  of  003(1/3)  at  4.3  pm  is  quite 
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Figure  30.  Time-Dependent  Effects  at  an  Altitude  of  100  km 
Illustrated  by  the  8°  and  20°  CVFs. 

different.  The  20°  CVF  shows  an  increase  in  radiance  over  the  8°  CVF 
due  to  the  resonant  energy  transfer  of  CO2  with  vibrationally  excited  N2 
(see  Subsection  5.3). 

The  effect  of  altitude  on  the  electron  beam-induced  emission  is 
shown  in  Figs.  31  and  32  for  the  8°  and  20°  CVFs,  respectively.  These 
data  were  taken  during  the  downleg  portion  of  the  flight.  For  the 
prompt  emission  from  NO,  quenching  not  only  becomes  increasingly 
important  at  longer  times,  but  also  becomes  more  important  as  altitude 
decreases.  However,  for  NO,  the  most  dominant  quenching  process  appears 
to  be  the  quenching  of  the  precursor  for  NO,  N(2d).  The  effects  of 
quenching  the  high  vibrational  states  of  NO  are  also  shown  at  long 
times.  For  CO2,  the  quenching  of  C02(^'3)  and/or  vibrationally  excited 
N2  is  also  evident.  On  the  other  hand,  CO(Av=l)  is  not  effectively 
quenched  and  shows  the  effect  of  a  relatively  long  vibrational  lifetime. 
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gure  31.  Altitude  Dependence  of  the  8°  CVF  IR  Data. 


Wovelength 

Figure  32.  Altitude  Dependence  of  the  20°  CVF  IB  Data. 
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5 . 3  EXCEDE  Kinetic  Model 


The  most  important  IR  emitters  in  an  electron-!  jed  atmosphere 
include  C02»  CO,  and  NO.  Emissions  from  O3  are  also  of  potential 
systems  interest.  However,  much  less  is  known  about  O3  emission  during 
disturbed  conditions,  and  therefore  it  is  much  harder  to  treat  its 
emission  in  detail.  We  have  developed  a  general  kinetic  model  (EXCEOE) 
to  describe  the  chemical  and  collisional  processes  occurring  during  the 
electron  dosing  of  the  atmosphere.  This  code  is  based  on  the  CHEMKIN 
code  developed  at  Sandia  Laboratories. The  computer  code  uses  a 
finite  rate  approach  in  modeling  the  generation  of  secondary  electrons 
and  the  resultant  molecular  emission  from  collisions  as  well  as  from 
chemiluminescent  reactions.  In  particular,  the  model: 

•  Solves  the  time  and  energy  dependent  rate  equations  for 
arbitrary  dose  rates  and  times  to  obtain  species  number 
densities, 

•  Treats  the  chemical  and  collisional  processes  for  each 
vibrational  state  of  interest, 

•  Assumes  a  simplified  histogram  treatment  of  the  secondary 
electron  spectra,  and 

•  Uses  an  optically  thin  line-by-line  treatment  of  molecular 
emission  with  the  rotational  temperature  equal  to  the 
eunbient  temperature  to  obtain  volumetric  emissions. 

Since  the  IR  emission  resulting  from  direct  electron  impact 
excitation,  as  well  as  recombination  reactions,  depends  on  the  low 
energy  secondary  electron  spectra,  it  is  important  that  the  histogram 
treatment  retain  the  major  features  of  the  secondary  electron  spectra. 
Secondary  electrons  are  formed  via  collisional  ionization  of  M(=  N2 ,  03, 
or  O) 

ep  +  M  -»  fip  +  products  .  (30) 

Prom  an  assumed  initial  secondary  electron  energy  distribution, 
secondary  electrons  undergo  further  collisions  with  major  species 
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(31) 


•b  +  M  -»  +  M* 

and 

Bg  +  M  -*  e^'  +  M'*’  (32) 

which  results  in  the  relaxation  or  cascading  of  the  secondary  electron 
spectra. 

5.3.1  Molecular  Radiators 

The  three  major  IR  radiators  presently  included  in  the  EXCEDE  code 
are  NO,  CO2,  and  CO.  Descriptions  of  the  kinetics  for  NO  and  CO2  are 
given  in  Subsection  3.1.2. 

CO  emission  in  the  4.7  ptn  region  is  important  in  the  atmospheric 
window  to  the  red  of  CO2  *^3-  The  primary  source  of  excitation  is  from 
secondary  electron  excitation.  Vibrationally  excited  CO  is  formed 

through  shape  resonance  excitation  by  secondary  electrons  in  the 
spectral  regions  of  interest.  Thus,  the  CO  signature  is  calculated  by 
the  same  technique  as  used  to  calculate  CO2  ‘^2'  primary  reaction 

set  used  in  these  calculations  is  as  follows: 

eg  +  CO  -♦  e^  +  CO(i')  (33) 

CO(i/)  -*  CO  ♦  hi/  .  (34) 

The  CO  vibrational  state  populations  predicted  by  Reactions  (33)  and 
(34)  are  in  excellent  agreement  with  those  deduced  from  analysis  of 
earlier  EXCEDE  flight  data. 

5.3.2  .Model  Predictions 

The  chemical  kinetics  mechanisms  detailed  above  have  been 

incorporated  into  the  EXCEDE  code  data  base.  Preliminary  calculations 
have  been  performed  as  a  function  of  dose  time  and  altitude  in  order  to 
compare  with  the  EXCi.  ^ J  III  CVF  data.  It  should  be  emphasized  that  the 
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model  calculations  are  preliminary  in  two  ways:  (1)  the  dose  rates  used 
are  obtained  from  the  ARCTIC^l  model,  as  the  beam  dose  rate  profiles 
from  the  nJ  first  positive  319.4  nm  emission  are  not  yet  available  from 
the  data  and  (2)  the  natural  background  has  not  been  included  in  the 
calculations.  Therefore,  only  a  qualitative  comparison  with  the  data 
should  be  made. 

The  time  dependence  of  the  spectral  emission  for  NO(Av=l,2), 
002(1^3),  and  CO(Av=l)  at  100  km  predicted  by  the  EXCEDE  kinetic  model  is 
shown  in  Fig.  33.  Recall  that  the  0°,  8°,  and  20°  CVFs  correspond 

approximately  to  0,  0.2,  and  0.5  sec  after  beam  irradiation  of  the 

atmosphere.  The  trends  in  the  model  prediction  should  be  compared  with 
the  data  presented  in  Fig.  30.  The  qualitative  agreement  between  the 
model  and  the  data  is  very  good.  For  No,  the  effect  of  quenching  as  a 
function  of  time  is  clearly  seen.  The  importance  of  both  direct 
electron  excitation  of  CO2  and  vibrational-to-vibrational  energy 
transfer  of  vibrationally  excited  N2  with  CO2  to  form  C02(*'3)  is  also 
evident.  From  the  calculations,  the  emission  at  0°  corresponds  to 
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Figure  33.  Model  Prediction  of  the  Time-Dependence  of  the  0°,  8°, 
and  20°  CVFs  at  an  Altitude  of  100  km. 
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direct  electron  excitation  (35),  while  8°  and  20°  correspond  to 
vibrational  energy  transfer  with  Nj  (Eqs.  19  and  20)  .  In  complete 
agreement  with  the  data,  the  model  shows  that  C02(4.3  /jm)  emission 
increases  with  time. 

The  calculated  altitude  dependence  of  the  electron  beam-induced 
emission  is  shown  in  Figs.  34  and  35  for  the  8°  and  20°  CVFs, 
respectively.  These  calculations  should  be  compared  with  the  data  shown 
in  Figs.  31  and  32.  Although  the  trends  illustrated  by  the  model 
calculations  are  consistent  with  the  data,  there  are  two  discrepancies 
worth  mentioning.  First,  the  model  predicts  that  the  NO  emission  pea)«s 
at  100  )un,  while  the  data  indicate  that  NO  emission  is  greatest  at  95 
km.  This  observation  suggests  that  the  quenching  rate  of  N(2d)  by  O 
(Eq.  15)  may  be  too  fast.  Indeed,  there  is  some  disagreement  between 
laboratory  measurements  of  the  rate  constant  and  the  rate  constant 
deduced  from  atmospheric  observations.  Secondly,  the  data  show  that  CO 


emission  decays 

at  a 

rate 

much 

slower 

than 

expected 

based  on  the 

vibrational  lifetime. 

Since 

the 

kinetic 

model 

for  CO 

only  includes 

direct  electron 

excitation. 

the 

calculated  CO 

emission 

is  extremely 

small  for  the  8°  and  20°  CVFs  (also  see  Fig.  33).  This  observation 
suggests  that  another  mechanism  for  CO(Av=l)  emission  may  exist. 

5.4  Summary 

Excellent  electron  beam-off  and  beam-on  infrared  data  have  been 
obtained  in  the  altitude  range  of  85  to  115  km.  The  use  of  three  CVFs 
offset  from  each  other  allows  the  measurement  of  time-dependent  effects 
in  an  electron-dosed  atmosphere.  Significant  enhancement  of  the 
quiescent  atmosphere  has  been  observed  for  NO  (2.7  and  5.3  pm),  CO2 
(4.3  pm),  and  CO  (4.7  pm).  Inconclusive  evidence  has  been  found  for  NO'*’ 
(4.3  pm),  O3  (9.6  pm),  and  CO2  (15  pm)  enhancement.  However,  it  should 
be  stressed  that  the  lack  of  significant  enhancement  for  certain  species 
merely  indicates  that  the  background  radiation  dominates  radiation  from 
p>ossible  local  electron-beam  excitation.  Comparison  of  the  data  with  a 
kinetic  model  developed  for  EXCEDE  indicate  good  qualitative  agreement 
for  both  altitude  and  time-dependent  processes. 
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6.  CONCEPTS  FOR  A  SMALL  ROCKET  EXPERIMENT 


Experimental  measurements  of  the  altitude  profiles  of  various 
ambient  species,  especially  atomic  oxygen,  will  help  resolve  questions 
and  uncertainties  in  current  models.  One  possible  measurement  concept 
which  merits  strong  consideration  is  the  utilization  of  small  sounding 
rockets  to  directly  measure  species  profiles.  This  concept  has  an 
advantage  in  that  launch  costs  would  be  substantially  reduced  while 
still  having  the  capability  to  perform  multiple  experiments.  Reduced 
costs  will  permit  multiple  launches  to  provide  information  about  species 
variability.  The  main  constraint  is  that  the  measurement  system  must  be 
small  enougn  to  fit  into  a  small  rocket  and  still  be  able  to  garner 
useful  data.  The  feasibility  of  such  measurements,  including  possible 
instrument  configurations,  were  studied  by  Telic  Optics,  Inc. hss, 
Inc. Telic  Optica  performed  generic  trade-off  and  stray-light 
analyses  for  an  optical  telescope,  and  HSS  investigated  the  feasibility 
of  packaging  a  high  resolution  spectrometer  into  a  small  volume.  This 
section  summarizes  some  key  results. 

6 . 1  Stray  Light  Analysis 

Local  IR  atmospheric  emissions  at  high  altitudes  are  several  orders 
of  magnitude  less  than  the  upweiling  radiation  from  the  lower  atmosphere 
and  earth  surface. 54  xhe  small  telescope  apertures  due  to  a  small 
diameter  vehicle  can  lead  to  severe  problems  with  stray  light 
contaminating  the  desired  observations.  Several  candidate  optical 
systems  for  a  small-rocket  measurement  system  were  evaluated  for 
stray-light  effects.  Small  aperture-to-wavelength  ratios  generic  to 
these  systems  cause  diffracted  stray  light  to  be  a  major  contributor  to 
system  stray  light,  even  though  such  systems  usually  have  both  a  field 
stop  and  a  Lyot  stop. 
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The  results  of  the  optical  trade-off  study  can  be  summarized  in 
graphs  of  the  Point  Source  Transmittance  (PST)  and  Stray  Light 
Equivalent  Spectral  Radiance  (SLESR).  The  general  conclusions  from  the 
Teiic  Optics  study  are:^^ 

1)  The  stray  light  is  dominated  by  diffracted  stray  light, 
even  at  shorter  IR  wavelengths  (3-4  pm); 

2)  The  L/D  ratio  of  the  baffle  must  be  about  12  or  greater 
to  achieve  the  needed  off-axis  performance,  that  is,  the 
aperture  stop  and  primary  mirror  must  be  shaded  from 
earthshine  at  all  look  angles;  and 

3)  Performance  can  always  be  improved  by  making  the  baffle 
longer  or  the  aperture  stop  larger. 


6.1.1  Prototype  Telescope  Baffle 

Stray  light  suppression  in  a  small-aperture  system  can  be 

dramatically  improved  by  making  a  baffle  that  is  oversize  in  diameter 

compared  to  the  aperture  atop  of  the  system.  When  the  off-axis  angle  to 

the  earth  is  large  enough  so  that  the  stop  is  not  illuminated,  any  stray 

light  must  make  four  diffractions  before  impinging  the  FOV.  First, 

light  is  diffracted  at  the  baffle  opening,  then  at  the  aperture  stop  (at 

the  primary  mirror),  then  at  the  field  stop,  and  finally  at  the  Lyot 

stop.  At  each  diffraction  the  irradiance  is  reduced  by  a  factor 

proportional  to  approximately  where  6  is  the  angle  by  which  the 

6  ^ 

light  diffracts  at  that  particular  aperture.  Oversizing  the  diameter  of 
the  baffle  opening  can  easily  increase  the  first  diffraction  angle  by  a 
factor  of  four  (that  is,  from  1°  to  4°),  thereby  reducing  the  multiple 
diffraction  contribution  to  stray  light  by  a  factor  of  64.  This  reduces 
the  diffracted  stray  light  to  a  level  comparable  to  the  baffle 
(multiple-scatter)  contribution  to  stray  light. 

The  oversize  baffle  concept  was  used  to  design  a  prototype  baffle 
system,  which  is  shown  in  Fig.  36.  This  was,  in  turn,  used  to  study  the 
fore-optics  of  a  small-aperture  spectrometer.  Parameters  used  for  the 
foreoptics  include  a  1:1  relay  consisting  of  confocal  paraboloids,  but 
with  the  aperture  stop  at  the  first  paraboloid,  giving  an  L/D  ratio  for 
the  stop  equal  to  the  L/D  ratio  for  the  primary  mirror.  In  comparison 
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Figure  36.  Small  Aperture  (One  Inch)  Foreoptics  Design. 

to  a  large  (four-inch)  aperture  system,  this  small-aperture  system  gives 
equal  performance  at  large  earth  angles  (greater  than  20°),  and  superior 
performance  at  small  earth  angles  (10  to  20°).  The  baffle  opening  is 
oversized  so  that  the  edge  of  the  baffle  opening  is  3°  away  from  the 
near  edge  of  the  primary  mirror. 

The  baffle  vanes  are  closely  spaced  so  that  no  first  diffuse 
reflections  strike  the  primary  mirror,  and  furthermore,  no  first  diffuse 
reflections  from  the  outer  baffle  vane  cavity  can  strike  the  unvaned 
portion  of  the  baffle  tube. 

6.1.2  Prototype  Optical  System 

One  feasible  optical  system  has  a  collecting  aperture  of  only  1" 
which  allows  for  the  use  of  IX  confocal  paraboloids.  A  wide  FOV  can  be 
obtained  by  correcting  the  unit-magnification  confocal  paraboloids  for 
all  aberrations  except  field  curvature. 

The  imaging  capability  of  the  telescope  is  felt  to  be  of  secondary 
importance  since  the  image  quality  needs  only  to  be  sufficient  to  gather 
light  from  within  the  FOV  and  focus  it  on  the  spectrometer  slit.  Of 
greater  importance  is  the  control  of  aberrations  between  the  aperture 
stop  and  Lyot  stop.  Diffraction  is  the  dominant  source  of  stray  light 
in  a  small-aperture  system,  and  it  is  important  that  rays  diffracted  at 
the  aperture  stop  do  not  leak  through  the  Lyot  stop  because  of  pupil 
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aberrations.  The  aperture  atop  is  the  edge  of  the  primary  mirror.  That 
edge  is  a  circle  tilted  with  respect  to  the  LOS.  Its  image  is  also 
tilted. 

The  aperture-stop  image  is  also  tilted  because  of  the  telescope 
field  curvature.  Tilting  the  Lyot  stop  increases  its  effectiveness. 
Figure  36  shows  a  schematic  of  this  proposed  optical  design,  and  Table  3 
lists  relevant  optical  parameters. 

The  PST  has  been  calculated  for  the  fjffle  and  foreoptics  described 
above.  The  PST  for  a  9.2-10.0  pm  bandpass  is  shown  in  Fig.  37.  The 
oversize  baffle  greatly  reduces  the  diffracted  component  of  the  PST  to 
the  point  where  it  is  small  compared  to  the  baffle  scatter  component  for 
angles  outside  of  the  direct  shading  angle. 

The  Stray  Light  Equivalent  Spectral  Radiance  (SLESR)  is  the  amount 
of  apparent  spectral  radiance  that  the  system  sees  due  to  the  presence 
of  stray  light  in  the  system.  It  can  be  calculated  from  the  PST  by 
integrating  the  spectral  radiance  over  the  extended  source  (the  earth) 
weighted  by  the  PST,  normalized  to  one  steradian.  The  reduced  PST 
results  in  SLESR  values  which  meet  a  nominal  sensitivity  goal  of 
10“12  w/cm^/sr/pm.  Predicted  values  of  SLESR  are  shown  in  Fig.  38  for 
this  baseline  25  mm  aperture  system. 

Table  3.  IX  Telescope  Optical  Prescription. 


(■■■ 


I  Element 

Radius 

Thickness 

Material 

Clear 

Aperture 

Diameter 

{  Baffle 

1  Opening 

Inf inite 

300  mm 

Air 

70.0  X  88.0 

1 

!  Primary 

i 

-127.0  cc 

-63.5 

Ref  1 

25.4 

1 

1  Field  Stop 

Inf inite 

-63.5 

Air 

2.382  X  7 . 153 

j  Secondary 

127.0  cc 

127.0 

Ref  1 

30.2  X  39.7 

!  Lyot  Stop 

Infinite 

- 

Air 

22.9 
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Point  Source  Transmittance 

j  3,  25mm  ap,  300mm  baf.,  4deg  HFOV 


-a 

o 

- 

^<0 

_o 

Off-Axis  Angle  (degrees) 


6. 2  Prototype  Instrumentation  Package 


HSS,  Inc.  developed  and  studied  a  conceptual  design  for  a  small- 
rocket  IR  grating  spectrometer.^^  Sensitivity  is  an  important  issue 
because  atmospheric  emissions  are  weak  in  the  window  regions. 
Requirements  for  a  grating  spectrometer  are  driven  by  constraints 
imposed  by  the  small  size  of  the  rockets.  The  advent  of  solid  state 
photomultiplier  (SSPM)  detectors  now  makes  it  possible  to  design  compact 
grating  spectrometers  that  can  achieve  the  needed  spectral  resolution 
and  NESR.  Further  cost  reductions  would  be  achieved  by  developing  a 
single  reusable  spectrometer  that  spans  the  full  4.5  to  22  /jm  range.  On 
any  given  flight  this  spectrometer  would  measure  the  atmospheric 
emissions  in  only  a  small  spectral  region  of  1-2  pm;  grating  and 
detector/slit  parameters  would  be  changed  for  each  flight  to  achieve 
optimum  performance  for  the  desired  spectral  region. 

Spectrometer  performance  goals  are  listed  in  Table  4.  Also  shown 
are  characteristics  of  a  spectrometer  that  is  likely  to  meet  these 
goals.  A  spectral  resolution  of  2  cm”^  is  assumed,  although  resolutions 
from  1  to  4  cm“^  are  undoubtedly  acceptable  depending  upon  tradeoffs 
necessary  to  achieve  other  objectives.  The  threshold  sensitivity  is 
equivalent  to  1/30  of  the  Zodiacal  background  radiance  (ZLB),  which  is 
about  3  X  10"^^  w/cm^/sr/pm  over  the  4-20  pm  region.  Requiring  the 
dynamic  range  of  the  spectrometer  to  be  1000:1  is  related  more  to 
characteristics  of  the  SSPM  detector  than  to  the  the  spectrometer 
design . 

6.2.1  Design  Constraints 

The  spectrometer  will  be  flown  on  a  rocket  having  a  diameter  of 
from  9  to  14  inches.  Using  a  SSPM  detector  demands  that  the 
spectrometer  be  housed  in  a  cryogenic  vacuum  dewar  with  an  outside 
diameter  of  7  to  12  inches,  depending  upon  the  rocket  diameter.  The 
dewar  constrains  the  outside  diameter  of  the  spectrometer  to  be  5  to  10 
inches,  depending  upon  the  size  of  the  rocket.  Total  allowable  payload 
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Table  4.  Design  Goals  for  a  Small  Rocket  Spectrometer . ^3 


PERFORMANCE 

operating  Range . 4.5  to  22  pm 

Spectral  Resolution  .  2  cm'l 

NESR  .  w/ (cm^ . sr -pm) 

S/N  Ratio . 30:1 

(Zodiacal  Light  Background) 

Dynamic  Range  .  103:1 

CHARACTERISTICS 

Type:  LHe-Cooled,  Prefiltered  Grating  Spectrometer 

Scan  Time . 2  sec 

Field-of-View  .  2  degrees 

Relative  Aperture  .  f/2.5 

Throughput . 6  x  10“3  cm^  sr 

Detector: 


Type . SSPM 

No.  Channels . 1 


Objective  Optics: 

High  Resolution  Off-Axis 
Telescope 
5  cm2 

High  Efficiency  Design; 
1/D  >  12 

Size . Compatible  with  a  9-14 

inch  diameter  rocket 

Weight  (Including  Cryostat)  .  100  to  200  lbs. 


Type  . 

Aperture 

Baffle 


will  be  in  the  vicinity  of  500  to  600  pounds  with  100  to  200  pounds 
allocated  to  the  spectrometer  and  cryogenic  dewar;  the  remainder  is 
allocated  to  payload  support  systems. 

The  LOS  of  the  spectrometer  will  be  oriented  by  the  payload  support 
system  so  that  its  FOV  lies  slightly  above  the  earth's  limb.  A 
telescope  and  baffle  system  is  required  to  sharply  define  the  FOV  and 
reject  off-axis  radiation  from  the  earth's  limb  and  lower  atmosphere. 
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The  telescope's  objective  will  be  about  one  inch  in  diameter  and  the  FOV 
will  be  about  2  with  a  throughput  of  6  x  10~3  cm^  sr.  Efficient 
coupling  of  the  telescope  to  the  spectrometer  is  required  to  maximize 
sensitivity. 


6.2.2  Diffraction  Grating  Selection 

The  requirement  for  high  linear  dispersion  in  a  spectrometer  is 
customarily  met  by  use  of  long  focal  length  optics  or  high  angular 
dispersion,  or  a  combination  of  both  approaches.  Long  focal  length 
optics  imply  large  instruments,  an  approach  necessarily  ruled  out  in  the 
present  situation  where  a  small  compact  spectrometer  is  required. 

The  alternative  approach,  high  angular  dispersion,  requires  working 
at  large  diffraction  angles.  At  large  diffraction  angles,  however,  wide 
spectral  ranges  are  not  possible  with  a  single  grating  spacing  unless  a 
coarse  ruling  is  emp  eyed.  This,  in  turn,  leads  to  large  spectral  order 
numbers,  where  many  orders  are  required  to  cover  a  given  spectral  range. 
At  extreme  diffraction  angles,  the  physical  length  of  the  grating 
becomes  several  times  the  aperture  dimension  so  that  for  large  apertures 
and  fast  f/numbers,  the  grating  size  becomes  excessive. 

The  HSS  study^3  concentrated  on  the  use  of  coarse  dirfraction 
gratings  and  moderate  diffraction  angles  to  avoid  an  excessive  grating 
size.  This  approach  necessarily  led  to  the  use  of  multiple  orders  to 
cover  the  spectral  range  4.4  to  22  pm. 

Nominal  grating  and  instrument  parameters  which  are  consistent  with 
all  other  aspects  of  the  spectrometer  design  are: 


Grating  Constant: 
Diffraction  Angle: 

Focal  Length: 

Detector  Element  Width: 


a  =  50  pm 
P  =  20° 
f  =  200  mm 
w'  =  0.150 


6.2.3  Spectrometer  Design 

Because  of  the  large  spectral  range  (4.4  to  22  pm)  to  be  covered  by 
the  spectrometer  and  the  lack  of  suitable  refractive  optical  materials 
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for  Buch  a  range,  reflective  rather  than  transmission  optics  are 
indicated.  Four  basic  spectrometer  configurations  were  considered:  the 
reflective  (Littrow),  the  Ebert,  the  Conventional  Czerny-Turner  (CCT) 
and  the  Crossed  Czerny-Turner  (XCT).  These  configurations  were  selected 
for  examination  because  of  their  known  qualities  of  simplicity, 
compactness  and  fast  relative  apertures. 

Primary  consideration  was  given  to  throughput  and  image  quality. 
The  Littrow  configuration  was  rejected  because  of  lack  of  first-hand 
experience  with  it  in  the  reflective  version.  (It  appears  to  be  more 
adaptable  to  refractive  rather  than  reflective  optics.)  The  Ebert 
configuration  was  rejected  in  favor  of  the  Czerny-Turner  configurations 
because  the  latter  configurations  offer  larger  numerical  apertures  with 
good  image  quality.  Finally,  the  XCT  configuration  was  chosen  over  the 
CCT  configuration  because  experience  has  shown  that  for  a  given 
numerical  aperture  the  XCT  has  smaller  off-axis  angles  than  the  CCT  and, 
hence,  better  resolution.  The  XCT  also  has  the  advantage  over  the  Ebert 
and  CCT  in  that  it  is  less  susceptible  to  the  problem  of  re-entry 
spectra. 

The  basic  shape  of  the  spectrometer  will  have  the  form  of  a 
cylinder  oriented  with  the  axis  of  the  cylinder  parallel  to  the 
longitudinal  axis  of  the  carrier  rocket.  As  described  earlier,  the 
radial  dimension  of  the  spectrometer  is  critical  in  that  it  is  dictated 
by  the  size  of  the  cryogenic  dewar.  The  diameter  of  the  dewar  is 
determined  in  turn  by  the  diameter  of  the  carrier  rocket.  The  possible 
range  of  9  to  14  inches  for  the  diameter  of  the  rocket  translates  into  a 
possible  range  of  5  to  10  inches  for  the  radial  dimension  of  the 
spectrometer.  There  is  no  critical  constraint,  within  reason,  on  the 
length  of  the  spectrometer. 

There  are  two  alternative  applications  of  the  basic  spectrometer 
design,  one  as  a  scanning  monochromator,  the  second  as  a  spectrograph. 
The  spectrograph  version  has  the  potential  of  stationary  optics  and  a 
multiplex-advantage.  The  scanning  monochromator  has  the  potential 
advantages  o  greater  throughput  and  superior  spectral  resolution.  Both 
versions  of  the  spectrometer  were  examined  at  some  length. 
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6.2.4  SSPM  Detectors 


Critical  to  the  performance  of  the  spectrometer  are  the 
capabilities  and  limitations  of  the  solid  state  photomultiplier  (SSPM) 
detectors  intended  for  use  with  the  spectrometer.  Information  of  this 
nature  was  supplied  by  technical  personnel  of  AFGL  and  the  SSPM 
development  team  at  the  Roc)cwell  International  Science  Center. 

SSPM's  are  small  area  (pixel  size)  photon  counting  detectors.  They 
are  a  relatively  new  development  in  detection  techniques  and  one  may 
expect  dramatic  advances  in  their  capabilities  in  the  coming  years.  The 
detectors  wor)c  well  at  cryogenic  temperatures,  7  to  10°K,  but  apparently 
do  not  work  well  at  the  temperature  of  liquid  helium  (4.20K).  They  have 
response  capability  extending  from  the  visible  to  approximately  JO  pm, 
with  a  quantum  efficiency  which  depends  both  on  photon  wavelength  and  on 
the  type  of  illumination,  front,  back  or  edge. 

For  a  small  rocket  spectrometer,  the  signal-to-noise  ratio  (S/N)  is 
perhaps  the  most  important  indicator  of  spectrometer  performance.  The 
noise  analysis  considered  two  primary  sources:  (1)  the  dark  count  of 
the  SSPM  and  (2)  photon  noise  associated  with  the  Zodiacal  Light.  A 
comparison  of  noise  predictions  for  two  versions  of  the  spectrometer^^ 
showed  that  the  spectrograph  has  a  average  factor  of  four  to  six 
advantage  in  S/N  over  the  monochromator,  che  exact  value  depending  on 
bandpass  and  specific  design  parameters. 
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7. 


IMPULSE  APPROXIMATION  FOR  ATOM-MOLECULE  COLLISIONS 


The  GL  Field-Widened  Interferometer  (FWI)  program^^  observed 
emissions  from  excited  states  of  NO  in  an  aurorally  disturbed 
atmosphere.  R-branch  band  heads  for  the  3-2,  2-1,  and  1-0 
vibrational-rotational  (VR)  transitions  were  observed  at  1960,  1990,  and 
2020  cm“l,  resp>ectively .  The  appearance  of  these  band  heads  leads  to 
the  identification  of  these  emissions  as  high  rotational  states. 
Modeling  studies  which  used  a  two-temperature  model  for  the  NO 
rotational  states  indicate  that  the  band  heads  do  not  become  observaDle 
until  the  high-temperature  component  reaches  2500  K.  Thus  the  NO  is  in 
a  highly  NLTE  state.  One  hypothesis  for  the  formation  of  these  states 
is  that  they  result  from  VR  transfer  of  energy  during  the  collision  of 
ambient  species  with  excited  atoms  which  have  been  produced  via 
chemiluminescent  processes. This  section  discusses  the  Impulse 
Approximation  (lA)  for  atom-molecule  collision  processes  and  its 
extension  to  include  the  effects  of  VR  energy  transfers. 

In  the  lA  three-body  process  are  treated  as  a  pair  of  atom-atom 
collisions,  with  the  third  atom  in  each  case  remaining  a  spectator. 
This  formalism  was  originally  developed  by  Chew58  the  context  of 
neutron-deuteron  collisions.  It  has  also  been  applied  to  atomic  physics 
problems  by  Bogan,®®  by  Eckelt,  Korsch,  and  Philipp,®®”®®  and  by  Beard 
and  Micha.®^  The  scattering  amplitudes  for  the  two  processes  are 
determined  by  the  momentum  space  wave  functions  of  the  molecular  initial 
and  final  states  and  the  two-body  T-matrices.  The  differential  cross 
section  is  determined  by  summing  over  the  two  two-body  amplitudes.  The 
earlier  work®®”®^  dealt  with  hyperthermal  collisions  where  the  energy  of 
the  projectile  atom  is  quite  high  (of  the  order  of  an  electron  volt  or 
higher).  In  Sharma's  suggestion,®®  low  projectile  energy  collisions  may 
be  treated  in  the  lA  when  the  molecular  internal  energy  is  large.  Thus 
it  is  the  relative  velocity  of  the  projectile  with  respect  to  the 
molecular  atom  that  matters  in  the  two-body  impulsive  collision.  A 
large  relative  velocity  is  achieved  by  having  either  a  fast  projectile 
or  a  highly  excited  molecule.  With  appropriate  modifications,  the  lA 
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formalism  may  then  be  capable  of  representing  a  large  variety  of 
collision  processes  that  are  applicable  to  atmospheric  processes. 

Previously  the  lA  has  been  employed  by  making  a  further 
approximation  called  the  Peaking  Approximation  (PA).®^"®^  The  typical 
two-body  process  is  represented  in  the  lA  by 

II  =  J*  0£(^3)  ipi(^3)  d^i  ,  (35) 

where  spectator  atom  momentum,  <p^  and  (pf  the  momentum  space 
wave  functions  of  the  initial  and  final  states,  and  T  is  the  two-body  T- 
matrix.  ^3  is  the  molecular  relative  momentum,  and  ^3  is  the  two-body 
relative  momentum,  with  the  relation 


^1  =  *1*3  “  ti^3 

(36) 

a  =  1  -  ib. 

mo 

b  =  — = — 

m2  m3 

t 

(37) 

where  1)3  is  the  incoming  (projectile)  momentum.  The  momenta  are  defined 
in  the  center-of-maas  reference  frame  for  the  collision  and  are 
express  ible®4  in  terms  of  any  two  of  the  Jacobi  momenta  and  For 
simplicity  we  describe  here  the  homonuclear  case  (mj  =  m2 ) .  Integration 
of  Eq.  (35)  is  quite  involved  and  is  simplif  ied^2-64  the  PA  by 
approximating  T  with  its  (fixed)  value  at  ^3  =  where  the  product  of 
the  two  wave  functions  is  peaked.  Then  T  can  be  pulled  out  of  the 
integral,  and  the  remaining  molecular  wave  function  overlap  integral 
leads  to  a  form  factor  which  depends  on  the  properties  of  the  target. 

While  the  lA  can  be  applied  to  many  physical  situations,  the  PA  is 
rather  restrictive.  In  order  to  employ  the  lA  for  a  wider  class  of 
phenomena,  it  is  thus  essential  to  develop  a  formalism  which  does  not 
invoke  the  PA.  Initial  steps  toward  such  a  formalism  was  developed  in 
the  course  of  this  work.  It  requires  a  more  involved  angular  momentum 
coupling  algebra  and  also  more  intricate  radial  integrals.  The 
directional  correlation  between  the  momentum  transfer  of  the  impulse 
collision  and  the  momentum  transfer  appearing  in  the  target  form  factor 
was  not  recognized  in  the  PA.  The  new,  more  general  formalism  fully 
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recognizes  this  correlation  and  can  determine  when  the  PA  is  (or  is  not) 
a  good  approximation.  In  terms  of  Eq.  (37),  it  can  be  shown  that 

if  is  dominated  by  the  projectile  momentum,  ^3,  the  PA  will  be 

reasonable  since  T  does  not  vary  much  over  the  domain  of  integration  in 
Eq.  (35).  On  the  other  hand,  for  processes  dominated  by  the  molecular 
momentum  ^3,  the  PA  will  not  be  adequate,  and  the  more  general  approach 
is  necessary. 

The  new  formalism  is  based  on  a  representation  for  the  T-matrix 
which  recognizes  the  invariance  of  the  momentum  transfer  Aj)  =  ^2  “  ?3  “ 
-  ^3  in  evaluating  the  integral  in  Eq.  (35).  The  T-matrix  can  be 
expanded  in  a  complete  set  of  spherical  harmonics^^ 

^  YLM(q3)  TL„(q3;Ap,p3,0)  ,  (38) 

LM 

where  0  is  the  angle  of  scattering  between  initial  and  final  projectile 
momenta,  ^3  and  ^3.  For  low  projectile  energies  (ap3  <<  bq3),  ^3  is 

essentially  along  ^3,  and  a  simpler  expansion  results, 

T(%*^i;e)  =  ^  PL(q3*^P)  TL(q3»'^P»®)  •  (39) 

L 

Both  Eqs.  (38)  and  (39)  are  found  to  be  rapidly  convergent  in  L,  and 
thus  this  Momentum  Transfer  Representation  (MTR)  offers  significant 
advantages  over  the  standard  expansion. The  expansion  coefficients 
Tlm  i-*'  Eq.  (38)  were  evaluated  numerically  for  the  hard  core  potential 
by  employing  the  orthogonality  relations  of  the  spherical  harmonics  and 
taking  the  T-matrix  for  the  hard  core  potential  as  given  in  the  standard 
representation  by  Van  Leeuwen  and  Reiner^^  and  Beard  and  Micha,^'' 

T(^',^;e)  =  ^  ^  (2^  t£(q',q;e)  P£(co3  0)  .  (40) 

e 

The  expansion  coefficients  Tl  in  Eq.  (39),  for  the  low  projectile  energy 
limit,  were  evaluated  numerically  by  employing  the  orthogonality 
relations  of  the  Legendre  functions.  For  the  "post"  form  of  the 
T-matrix,  characterized  by  c  =  h2q^/2p3^3,  it  was  discovered  from  the 
numerical  evaluations  that  Tl  had  the  exact  factorization  property 
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TL(q.^p) 


jL(^P)  tL(q) 


(41) 


This  remarkable  result  offers  considerable  computational  advantage  in 
subsequent  calculations.  It  can  also  be  shown  that  t£^(q)  are  related  to 
the  forward  scattering  amplitude  and  its  derivatives. 

In  the  new  formalism,  the  integration  in  Eg.  (35)  is  accomplished 
without  making  any  approximations.  This  is  a  significant  step  in  the 
development  of  the  lA,  which  can  now  be  applied  to  a  wider  variety  of 
physical  situations  since  it  is  no  longer  necessary  to  invoke  the  PA. 
The  general  approach  is  based  on  taking  ^3  as  the  variable  of 

integration  in  Eq.  (35).  T  is  expressed  in  terms  of  ^3  through 

Eq.  (38).  The  final  state  wave  function,  which  depends  on  ^3  =  ^3  - 
is  represented  through  an  inverse  Fourier  transform  in  terms  of  a 
coordinate  space  integral.  The  ^3  dependence  is  now  explicit,  and  the 
angular  and  radial  integrals  in  ^3  space  can  be  carried  out.  The 

angular  integral  gives  various  Clebsch~Gordan  (CG)  coef f icients. 
Computer  programs  were  developed  to  evaluate  the  radial  integrals  and 
the  various  CG  and  Racah  coefficients.  The  hyperthermal  limit  (ap3  » 
bq3)  is  compared  to  previous  treatments^O"^'*  for  which  the  PA  emerges  as 
limiting  cases  of  this  formalism.  The  precise  errors  involved  in 
invoking  the  PA  for  a  given  set  of  physical  parameters  can  be  determined 
with  this  formalism.  As  an  application,  this  approach  was  evaluated  by 
looking  at  the  cross  sections  for  N2  excitations  resulting  from 

collisions  with  Li'*’.^^  Figure  39  presents  an  illustrative  comparison  of 
the  calculated  differential  cross  sections  for  the  N2,  V  =  2  and  3, 
excited  scates  and  shows  that  the  PA  significantly  overestimates  the 
exact  lA  cross  section  for  this  excitation  processes.  Thus  the  PA  is 
inadequate  in  the  low  energy  limit,  and  this  new  formalism  offers  a 
uniform  method  for  any  projectile  and  molecular  energies. 

Possible  extensions  for  the  next  higher  level  of  approximation  are 
the  double  collision  terms  in  the  expansion  of  the  scattering  matrix. 
Three  feasible  approaches  that  merit  further  study  are; 
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Figure  39.  Differential  Cross  Sections  (ctn^)  for  the  Exact  and  PA 
Calculations  vs.  Energy  Loss  (eV)  for  the  Collision 
Process  Li"*"  +  N2(v=0)  -+  Li"*"  +  N2(v=2,3)  at  a 

Center-of-Mass  Scattering  Angle  of  30°  and  an  Initial 
Translational  Energy  of  1  eV.  The  Peaking  Calculations 
(o  and  A)  give  Larger  Cross  Sections  than  the  Exact 
Calculations  (Solid  Symbols).  (From  Ref.  58.) 


1)  Calculate  each  process  (p-*p'',  p"-»p’ )  in  its  momentum  transfer 
representation  (z  along  q'=p"-p  and  q"=p'-p",  respectively).  This 
approach  allows  the  individual  collision  amplitudes  to  be  calculated 
with  minimal  computation.  The  integration  over  the  intermediate 
momentum  p",  and  the  summation  of  the  intermediate  molecular  states  is 
not  practical  in  this  approach  as  it  would  require  numerical 
integration. 

2)  Another  approach  is  to  take  z,  the  quantization  axis,  along  p, 
the  incoming  projectile  direction.  The  second  process  can  be  computed 
using  as  the  quantization  axis.  While  the  individual  amplitudes  can 
again  be  written  down  and  computed  (with  slightly  more  computational 
effort  than  in  case  1),  the  intermediate  integration  over  ]^"  is  still 
not  feasible,  analytically,  and  would  require  numerical  integration.  As 
in  case  1,  the  accuracy  of  the  final  result  would  be  questionable, 
unless  a  very  large  numerical  grid  is  employed. 
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3)  In  view  of  these  difficulties,  which  could  not  be  overcome  so 
far,  we  return  to  a  previous  strategy,  that  is  using  the  general 
expression  for  an  arbitrary  direction  of  quantization,  and  expressing 
the  single  collision  amplitude  in  terms  of  a  bilinear  harmonic  form  in 
terms  of  the  incoming  and  outgoing  directions.  Then  the  intermediate 
momentum  integration  over  the  angles  is  analytically  feasible.  This 
allows  one  to  express  the  double  collision  amplitude  as  an  integral  over 
the  magnitude  p"  and  the  summation  of  the  intermediate  molecular  states. 

In  terms  of  accuracy,  which  is  essential,  since  the  double 
collision  term  is  a  correction,  the  third  approach  may  be  the  best.  It 
is  computationally  extensive,  and  it  may  be  best  to  work  with  small  j 
and  j'  to  limit  the  computing  time.  Considering  the  same  processes  as 
in  the  cited  papers  (0,0-*a,6  or  0,2-*l,6,  etc.),^^'^^  this  would  allow 
one  to  determine  if  the  double  collision  term  can  extend  the  domain  of 
semi  detailed  balance  non-violation  to  smaller  scattering  angles. 
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8.  FINE-STRUCTURE  TRANSITIONS  IN  COLLISIONS  OF  OXYGEN  ATOMS 


A  detailed  formulation  of  the  scatterii*g  of  oxygen  atoms  in 
collisions  with  each  other  could  provide  a  basis  for  inferring  oxygen 
altitude  profiles  from  field  measurements  of  its  various  ground-state 
transitions.  As  part  of  this  formulation,  the  correspondence  between 
atomic  molecular  states  is  worked  out.  An  expression  is  obtained  f.'^r 
the  cross  section  for  transitions  amongst  the  fine-structure  states. 
The  molecular  potential  energy  curves  that  enter  into  the  cross  section 
calculations  are  assembled. 

8 . 1  Introduction 

The  interpretation  of  the  fine-structure  emissions  of  atomic  oxygen 
in  the  atmosphere  involves  a  determination  of  the  populations  of  the 
fine-structure  levels.  At  !ow  altitudes,  collisions  will  maintain  the 
levels  in  thermal  equilibrium,  but  at  some  critical  altitude  that 
depend'j  on  the  magnitude  of  the  collision  cross  sections  the  levels  will 
fall  out  of  thermal  equilibrium,  and  a  detailed  solution  of  the 
equations  of  statistical  equilibrium  is  required.  The  cross  sections 
are  uncertain.  The  only  calculation,  that  of  Allison  and  Burke, was 
based  on  a  first  order  perturbation  theory  of  the  interactions,  which 
cannot  be  accurate  at  the  low  velocities  at  which  collisions  take  place 
in  the  atmosphere.  Here  we  present  the  full  molecular  theory.  It 
involves  a  much  more  complicated  angular  momentum  coupling  process. 

8 . 2  Theoretical  Formulation 

The  total  Hamiltonian  may  be  written 


Ht 


Ha  +  Hb  + 


^  V  (  r  ig  ,  r  ^ b ) 
i-l 


+  Hfj 


(42) 
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where 


Ha 


i 


(43) 


is  the  electronic  Hamiltonian  of  one  of  the  atoms  a,  Hjg  is  the 
electronic  Hamiltonian  of  the  other  atom  b,  V(ria,rj)3)  is  the 
electron-electron  repulsion  between  the  electrons  in  a  and  b,  and  Hjj  is 
the  nuclear  Hamiltonian  containing  the  electrostatic  interactions 
between  the  electrons  of  each  atom  with  the  nuclei  of  the  other.  We 
write  the  eigenstates  of  H^  as  4>a( j ,  a»  •  •  •  ^N,  a )  H^  as 

djj ( rN+i ,  b '  •  •  •  •^2N, b)  •  ^e  designate  the  quantum  numbers  of  the  states  of 
a  by  Q  =  ( f a"'a®a‘*'sa)  where  iamg  are  the  orbital  angular  momentum  and 

projection  quantum  numbers  of  a,  and  similarly  we  define  =  ( ^^mj^sj^mg^ ) 
for  b . 

At  large  internuclear  distances  R,  the  molecular  eigenstates  are 
linear  combinations  of 

I  ■^a'^a^  1  ^b'^b^  I 

wheri;  S  and  Mg  are  respectively  the  total  spin  and  spin  projection 
quantum  numbers. 

It  is  convenient  to  introduce  the  angular  functions 


iLASMg> 


(-l)«b  ^  [mg  -a)  I  ^  an'a>  I  ^b"'b>  I  SMs>  (44) 

mamb 


which,  however,  are  not  eigenfunctions  of  the  total  angular  momentum 
L=ia+£b»  because  the  jtets  1  £ama> | J^b'^'b^  defined  on  different  centers. 

It  may  be  shown  by  following  procedures  used  by  Knipp,^^  Chang,®®  and 
fminskij  and  Ni)citin,®®  that  for  singlet  molecular  states  S=0,  the 
gerade  (g)  Born-Oppenheimer  states  separate  to  L=0  or  L=2  atomic  states, 
and  ungerade  (u)  states  separate  to  L=1  atomic  states.  For  triplet 
molecular  states  S=l,  the  g  states  separate  to  L=1  and  the  u  states  to 
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L=0  and  2  states,  and  the  u  states  to  L=l.  The  detailed  connections  are 
given  in  Table  5.  There  are  81  molecular  states  formed  in  the  approach 
of  two  oxygen  atoms. 

There  is  a  one-to-one  correspondence  between  the  atomic  and 
molecular  states  except  for  the  b^Sg  and  2^7:'^  states,  the  and  2^7:* 
states,  and  the  l^Sg  and  S^Tg  states  which  are  different  linear 
combinations  of  | L=0,M=0, SMa>  and  |l=2 ,M=0, SMa>  states  with  S=0,  1  and 
2,  respectively.  To  determine  the  correct  linear  combinations,  we  must 
consider  the  nature  of  the  interaction  at  large  internuclear  distances. 


Table 

5. 

Connections  between  Atomic 

and  Molecular 

States 

for  Two 

Oxygen  Atoms. 

NO.  OF 

D£)GEN  • 

INVERSION 

MOLECULAR 

STATES 

A 

S 

STATE 

ATOMIC  STATE 

2 

±2 

0 

g 

^Ag 

L=2  M=±2  S=0  Ms> 

2 

±1 

0 

g 

'”g 

L=2  M=±l  S=0  Ms> 

1  2 

±1 

0 

u 

lliu 

L=1  M=il  S=0  Ms> 

1 

0 

0 

u 

Iv 

L=1  M=0  S=0  Ms> 

2 

0 

0 

g 

bir^,  2i2^  { 

L=0  M=0  S=0  Ms>, 

jL=2  M=0  S=0  Ms>} 

6 

±2 

1 

U 

3Au 

L=2  M=±2  S=1  Ms> 

6 

±1 

1 

U 

L=2  M=±l  S=1  Ms> 

6 

±1 

1 

g 

3ng 

L=1  M=tl  S=1  Ms> 

3 

0 

1 

g 

L=1  M=0  S=1  Ms> 

6 

0 

1 

u 

a3j;^,  { 

L=0  M=0  S=1  Ms>, 

|l=2  M=0  S=1  Ms>> 

10 

±2 

2 

g 

^Ag 

L=2  M=J:2  S  =  2  Ms> 

10 

±1 

2 

g 

5ng 

L=2  M=tl  S=2  Ms> 

10 

±1 

2 

u 

^Hu 

L=1  M=±l  S=2  Ms> 

1  5 

0 

2 

u 

L=1  M=0  S=2  Ms> 

10 

0 

2 

g 

L=2  M=0  S=2  Ms>, 

|l=o  M=0  S=2  Ms>} 
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Expanding  the  interaction  potential  in  multiples 


OD 


Q=3 


where 


vq 


-  I 


Z 


Q-2 

Z  Z 

L=1  M 


X 


(45) 


(46) 


i  refers  to  electrons  on  atom  a  and  j  to  electrons  on  atom  b,  and  G(L,M) 
is  a  coefficient  given  by  Chang. For  ^P-^P  interactions  of  neutral 
atoms,  the  leading  term  is  the  quadrupole-quadrupole  interaction  Q=5. 
After  some  angular  momentum  algebra,  we  obtain  for  the  asymptotic  matrix 
element  of  the  interaction  potential 


<Ma'Mb-|V5lMaMb> 


ii 

n 


2  2 

R(Ma,Ma.  ;Mb,Mb'  ) 


(47) 


where  r^  is  the  expectation  value  for  a  single  p  electron  of  oxygen,  and 
R  is  a  sum  of  products  of  3-j  symbols 


R(Ma,Ma' ;Mb,Mb' )  =  ( -1 )”a+Ma • +Mb^Mb •  ^G(L=2,M) 

M 

X  ^  (-1)Mi--^M3-  (48) 


Ml, M3 
Ml- ,M3. 


f  1  2  1 

1  1  1 

1  1  1  1 

[  1 

2  1  1 

[-M3 '  M  Ml 

M3  M-  -Mbj 

M3’  M-  -Mb'J 

[-M3' 

-M  M3J 

The  matrix  element  of  V5  can  be  expressed  more  simply  by  introducing  a 
6-j  coefficient  in  the  form 
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<Ma.Mb.  lV5lMaMb>  =  ^  |{i  1  ijl^  ^  G(2,M)  (-l)”a'+Mb’  (49) 

M 

(-Ma  Ma-  -m]  (-Mb  Mb'  m) 

Now  define  states  of  angular  momentum  J=L+S  according  to 

|jMj>  =  VIjTI  2  |L=l,Ma>  lL=l,Mb>  .  (50) 

MaMb 

Then 

|V5|jMj>  =  ^  I  (i  i  l)l^  Z  g(2'M)  -V2J+1  -V2J-+1  (51) 

M 


X 


2  (wb  Ma  -Mj)  [ 


1 

Mb' 


1 

Ma- 


Ma.Ma' 

Mb, Mb' 


1 

Mb' 


Now 

g(2,M) 


Art  _ 4J _ 

5  (2-M)!(2+M)! 


and 


Hence 


<J*Mj  IVsl JMj> 


n 


1 

36 


<r^>2 

r5 


X(J'Mj. , JMj) 


(52) 


(53) 


(54) 


where 


M  Ma,Ma. 


f  1  1 

J  ] 

f  1 

1 

J'  ] 

[Mj-Ma  Ma 

-MjJ 

[-Ma'+Mj’ 

Ma- 

-Mj-  J 

(-Ma  Ma-  ”)  (wj-Ma  "Ma’W-  «) 


Numerical  evaluation  yields 


X(2222) 


X(2121)  = 


4  6 

--  ,  X(2020)  =  - 


All  other  matrix  elements  vanish  and  there  is  no  mixing  between  the 
|j=2,M=0>  and  |j=0,M=0>  states.  The  element  <0000>  vanishes,  and  <  020> 
is  positive.  So  |j=2,M=2>  is  associated  with  the  stat.,  and 
|j=2,M=0>  with  the  b^Zg  state. 

The  asymptotic  interaction  energies  of  the  molecular  states  have 
<r2>2 

the  form  C5 - - —  with  the  values  of  C5  listed  in  Table  6. 

We  now  summarize  the  angular  momentum  coupling  for  the  scattering 
problem.  Let  i  be  the  angular  momentum  of  the  relative  motion  of  the 
two  atoms.  We  couple  £  to  the  internal  angular  momentum  j=Ja+Jb  of  the 
atoms  to  form  j .  Introduce  scattering  wave  functions  with  relative 
momentum  k 


♦(k,R)  =  ^ 

JMj£j 


F[R;JMj£j: 


where 


|jMj£j>  =  (-l)j-^-Mj  -4Uri  ^  Y£„^(k)  |jM. 


>  .  (58) 
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Table  6.  Values  of  the  Parameter  C5 


MOLECULAR 

STATE 


ATOMIC  STATE 
M  S 


SINGLETS 


In. 


21Ec 

b^Er 


TRIPLETS 


2^2 


QUARTETS 


=Ac 


2 

2 

1 

1 

2 

0 

2 

2 

1 

1 

2 

0 

2 

2 

1 

1 

2 

0 


±2 

±1 

±1 

0 

0 

0 

±2 

tl 

±1 

0 

0 

0 

*2 

±1 

tl 

0 

0 

0 


0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 


csr  ! 


6 

-24 

0 

0 

36 

0 

6 

-24 

0 

0 

36 

0 

6 

-24 

0 

0 

36 

0 


Units  of 


25 
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Following  and  slightly  generalizing  the  procedures  used  for  example  by 
Roueff  and  Dalgarno,'^^  we  obtain  for  the  interaction  Hamiltonian  in  the 
molecular  description 


■  j'  |v|  jjeMjJ>  =  (-l)i  +  3‘-«-^ '-Mj-Mj*  V2£  +  l  V2je'+1 


Z  [J' S' -^]  (j  s -s)  Z  I  Z  Z 


LS  L’S’  MlMs  Ml- Mg. 


■La  Sa  Ja' 

Lb  Sb  Jb 

*  ‘ 

‘L  S  j  ' 

<L’ML.S’Ma 

•  |h 

•La'  Sg-  Jg-V 

Lb-  Sb'  Jb'l  (ml  Mg  -A]  (hl'  Mg.  -M  ) 
^  L'  S’  j  • 


Various  selection  rules  are  embedded  in  this  algebra  which  simplify 
the  sets  of  coupled  equations  that  have  to  be  solved  to  determine  the 
scattering  matrix.  We  are  at  the  state  of  deciding  upon  the  best 
representations  of  the  potential  energy  curves  which  enter  as  diagonal 
elements  of  the  interaction  matrix.  We  have  obtained  satisfactory 
representations  of  1^11^,  and  states  and  acceptable 
estimates  of  the  remainder.  These  final  results  need  to  be  tested  for 
their  sensitivity  to  the  adopted  potentials. 
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APPENDIX  A 


SHARC  CHEMICAL  KINETICS  SCHEME 


The  current  version  of  SHARC-2  contains  chemical  kinetics/react  ion 
mechanisms  for  ambient  infrared  radiance  from  CO2,  O3,  NO,  H2O,  and  CO. 
The  original  SHARC  radiation  (i.e.,  chemical  kinetics  mechanisms  and 
associated  rate  constants)  mechanisms  were  based  on  those  used  in  the 
HAIRM  code.^1  In  order  to  be  consistent  with  limited  SHARC  and  field 
data  comparisons,  as  well  as  recent  laboratory  measurements,  the 
mechanisms  and  associated  rate  constants  from  the  HAIRM  code  have  been 
extensively  revised. 

In  the  following  sections,  the  vibrational  excitation/deexcitation 
mechanisms  and  radiative  relaxation  processes  are  listed.  The  species  M 
represents  N2  and  O2,  while  M’  represents  N2,  O2,  and  O.  Also  included 
in  the  tables  (Tables  A1-A5)  are  the  rate  constants  associated  with  the 
kinetics  mechanisms,  as  currently  used  in  SHARC.  When  available, 
references  for  the  rate  constants  in  each  Table  are  also  provided.  Much 
of  the  work  on  rate  constant  measurements  prior  to  1974  (and  to  a  lesser 
extent  calculations)  has  been  reviewed  by  Taylor. 
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TABLE  Al.  REACTIONS  AND  RATE  CONSTANTS  FOR  CO2  MODEL. 


Reaction 


Rate  Constant*  Reference 


CO2 (01101) 

+  M  s  002(00001) 

+  M 

1. 50x10“22x3 

a 

CO2(01101) 

+  0  =  002(00001) 

+  0 

1.00xl0“13xl/2 

b 

002(10002) 

+  M  =  002(01101) 

+  M 

3.00x10-22x3 

c,d 

002(10002) 

+  0  s  002(01101) 

+  0 

2.00x10-12t1/2 

c,d 

002(02201) 

+  M  »  002(01101) 

+  M 

6.00x10-22x3 

c 

002(02201) 

+  0  =  002(01101) 

+  0 

4.00x10-13x1/2 

c 

OO2 ( 10001) 

+  M  =  002(01101) 

+  M 

3.00x10“22t3 

c,d 

OO2 ( 10001) 

+  0  a  002(01101) 

+  0 

2.00x10-13x1/2 

c,d 

002(02201) 

+  M  =  002(10002) 

+  M 

2.90x10-13 

e 

002(10001) 

+  M  a  002(02201) 

+  M 

5.80x10-13 

e 

OO2 ( 10001) 

+  M  s  002(10002) 

+  M 

1.70x10-13 

e 

002(11102) 

+  M  a  002(10002) 

+  M 

1.50x10-22x3 

c,d 

002(11102) 

+  0  a  002(10002) 

+  0 

1.00x10-13x1/2 

c,d 

002(11102) 

+  M  a  002(02201) 

+  M 

3.00x10-22x3 

c,d 

002(11102) 

+  0  a  002(02201) 

+  0 

2.00x10-13x1/2 

c,d 

002(11102) 

+  M  a  002(10001) 

+  M 

1.50x10-22x3 

c,d 

002(11102) 

+  0  a  002(10001) 

+  0 

1.00x10-13x1/2 

c,d 

002(03301) 

+  M  a  002(10002) 

+  M 

4.50x10-22x3 

c,  d 

002(03301) 

+  0  a  002(10002) 

+  0 

3.00x10-13x1/2 

c,d 

002(03301) 

+  M  a  002(02201) 

+  M 

9.00x10-22x3 

c,d 

002(03301) 

+  0  a  002(02201) 

+  0 

6.00x10-13x1/2 

c,d 

002(03301) 

+  M  a  002(10001) 

+  M 

4.50x10-22x3 

c,  d 

002(03301) 

+  0  a  002(10001) 

+  0 

3.00x10-13x1/2 

c,d 

002(11101) 

+  M  a  002(10002) 

+  M 

1.50x10-22x3 

c,d 

CO2 (11101 ) 

+  0  a  002(10002) 

+  0 

1.00x10-13x1/2 

c,d 

CO2 ( 11101) 

+  M  a  002(02201) 

+  M 

3.00x10-22x3 

c,d 

002(11101) 

+  0  a  002(02201) 

+  0 

2.00x10-13x1/2 

c,d 

002(11101) 

+  M  a  002(10001) 

+  M 

1 . 50x10-22x3 

c,d 

002(11101) 

+  0  a  002(10001) 

+  0 

1.00x10-13x1/2 

c,d 

TABLE  Al.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


CO2 (03301) 

+  M  a 

CO2(11101) 

+  M  = 

C02( 11101) 

+  M  = 

002(00011) 

+  M  5^ 

CO2 (00011) 

+  0  = 

CO2 (00011) 

+  M  a 

CO2 (00011) 

+  M  a 

CO2 (20003) 

+  M  = 

002(20003) 

+  0  s 

002(12202) 

+  M  s 

002(12202) 

+  0  = 

002(20002) 

+  M  s 

002(20002) 

+  0  » 

002(04401) 

+  M  = 

002(04401) 

+  0  - 

002(12201) 

+  M  a 

002(12201) 

+  0  « 

002(20001) 

+  H  ^ 

002(20001) 

+  0  - 

002(01111) 

+  M  = 

002(01111) 

+  0  s 

002(20003) 

+  M  s 

002(12202) 

+  M  s 

002(20002) 

+  M  » 

002(11102) 

+  M 

002(11102) 

+  M 

002(03301) 

+  M 

002(03301) 

+  M 

002(03301) 

+  0 

002(11102) 

+  M 

002(11101) 

+  M 

002(03301) 

+  M 

002(03301) 

+  0 

002(03301) 

+  M 

002(03301) 

+  0 

002(03301) 

+  M 

002(03301) 

+  0 

002(03301) 

+  M 

002(03301) 

+  0 

002(03301) 

+  M 

002(03301) 

+  0 

002(03301) 

+  M 

002(03301) 

+  0 

002(00011) 

+  M 

002(00011) 

+  0 

002(00011) 

+  M 

002(00011) 

+  M 

002(00011) 

+  M 

2.90x10“13  f 

1.70x10'13  f 

2.90x10~13  f 

1.00xl0“15  + 

5.16xlO-lle-76--75/Tl/3  ^ 

1.10x10-13^1/2  h 


1.00x10-13  + 

5. 16xl0-lle"’^^-'^^/1^'^^  d 
1.00x10-13  + 

5.16xl0-lle“'^®''^5/l‘^^^  d 


1.20x10-21t3  c,d 
8.00x10“13t1/2  c,d 
1.20x10-21t3  c,d 
8.00x10-13x1/2  c,d 
1.20x10-21x3  c,d 
8.00x10-13x1/2  c,d 
1.20x10-21x3  c,d 
8.00x10-13x1/2  c,d 
1.20x10-21x3  c,d 
8.00x10-13x1/2  c,d 
1.20x10-21x3  c,d 
8.00x10-13x1/2  c,d 
1.50x10-22x3  c,d 
1.00x10-13x1/2  c,d 
5.00x10-13  + 


2.58xl0-lle“‘^®-'^3/xl/3  ^ 

1.00x10-13  + 

5.16xl0-lle"'^®-'^^/1'^'^^  i 
1.00x10-13  + 

5.16xl0-lle-76-'^5/xl/3  . 


,-3 


TABLE  Al.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


CO2 (04401) 

+ 

M 

002(00011) 

+ 

M 

1.00x10-15  + 

5.16xl0-lle-'^6.75/Tl/3 

i 

CO2 (20001) 

+ 

M 

002(00011) 

+ 

M 

5.00x10-15  + 

2.58xlO-lle-'76.75/Tl/3 

i 

002(12201) 

+ 

M 

- 

002(00011) 

+ 

M 

1.00x10-15  + 

5.16xlO-lle-'^5.75/Tl/3 

i 

CO2 (04401) 

+ 

M 

002(20003) 

+ 

M 

2.90x10-13 

f 

CO2 (04401) 

M 

002(20002) 

+ 

M 

2.90x10-13 

f 

002(20001) 

+ 

M 

= 

002(04401) 

+ 

M 

5.80x10-13 

f 

002(04401) 

-f 

M 

= 

002(12202) 

+ 

M 

2.90x10-13 

f 

002(12201) 

M 

s 

002(04401) 

-f 

M 

2.90x10-13 

f 

002(12201) 

+ 

M 

002(12202) 

M 

1.70x10-13 

f 

002(20002) 

M 

» 

002(20003) 

+ 

M 

1.70x10-13 

f 

002(20001) 

+ 

M 

s 

002(20003) 

+ 

M 

1.70x10-13 

f 

002(20001) 

+ 

M 

a 

002(20002) 

+ 

M 

1.70x10-13 

f 

002(02211) 

-t- 

M 

a 

002(10012) 

+ 

M 

2.90x10-13 

f 

CO2 ( 10011) 

M 

a 

002(02211) 

+ 

M 

5.80x10-13 

f 

002(10011) 

+ 

M 

= 

002(10012) 

M 

1.70x10-13 

f 

002(03311) 

+ 

M 

002(11112) 

+ 

M 

2.90x10-13 

f 

002(11111) 

+ 

M 

002(03311) 

+ 

M 

2.90x10-13 

f 

002(11111) 

+ 

M 

002(11112) 

+ 

M 

1.70x10-13 

f 

002(04411) 

+ 

M 

= 

002(20013) 

+ 

M 

2.90x10-13 

f 

002(20012) 

M 

002(04411) 

+ 

M 

5.80x10-13 

f 

002(20011) 

+ 

M 

a 

002(04411 ) 

M 

5.80x10-13 

f 

002(20012) 

+ 

M 

002(20013) 

+ 

M 

1.70x10-13 

f 

002(20011) 

M 

002(20012) 

+ 

M 

1.70x10-13 

f 

002(20011) 

M 

a 

002(20013) 

+ 

M 

1.70x10-13 

f 

002(00011) 

+ 

N2(0)  =  002(00001)  + 

N2(1) 

8.91x10-12/t1/2 

g 

002(01111) 

+ 

N2(0)  -  002(01101)  + 

N2(1) 

8.91x10-12/t1/2 

j 

002(10012) 

N2(0)  •=  002(10002)  + 

N2(1) 

8.91x10-12/t1/2 

j 

TABLE  Al.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


CO2 (02211) 

+ 

N2(0)  =  002(02201) 

4 

N2( 

1) 

8.91x10-12/x1/2 

j 

C02( 10011) 

N2(0)  =  002(10001) 

4 

N2( 

1) 

8.91x10-12/t1/2 

j 

C02(11112) 

+ 

N2(0)  =  002(11102) 

4 

N2( 

1) 

8.91x10-12/t1/2 

j 

C02(03311) 

N2(0)  =  002(03301) 

4 

N2( 

1) 

8.91x10“12/t1/2 

j 

002(11111) 

+ 

N2(0)  =  002(11101) 

4 

N2{ 

1) 

8.91x10"12/t1/2 

j 

C02( 20013) 

+ 

N2(0)  -  002(20003) 

4 

N2( 

1) 

8.91x10"12/x1/2 

j 

CO2 (04411) 

+ 

N2(0)  =  002(04401) 

4 

N2( 

1) 

8.91x10~12/x1/2 

j 

C02( 20012) 

+ 

N2(0)  =  002(20002) 

4 

N2( 

1) 

8.91x10~12/x1/2 

j 

002(20011) 

+ 

N2(0)  c  002(20001) 

4 

N2( 

1) 

8.91x10"12/x1/2 

j 

002(01101) 

002(00001) 

4 

hi^ 

1.50 

k 

002(11102) 

002(00001) 

4 

hu 

0.001 

k 

002(11101) 

002(00001) 

4 

hi/ 

0.010 

k 

002(00011) 

002(00001) 

4 

hi/ 

432.95 

k 

002(10012) 

002(00001) 

4 

hi/ 

11.15 

k 

002(10011) 

002(00001) 

4 

hi/ 

17.92 

k 

002(20013) 

002(00001) 

4 

hi/ 

0.150 

k 

002(20012) 

-♦ 

002(00001) 

4 

hi/ 

0.718 

k 

002(20011) 

002(00001) 

4 

hi/ 

0.233 

k 

002(10002) 

002(01101) 

4 

hi/ 

1.15 

k 

002(02201) 

-4 

002(01101) 

4 

hi/ 

3.03 

k 

CO2 ( 10001) 

-4 

002(01101) 

4 

hi/ 

1.57 

k 

002(01111) 

-4 

002(01101) 

4 

hi/ 

423.35 

k 

002(11112) 

002(01101) 

4 

hi/ 

10.52 

k 

002(11111) 

002(01101) 

4 

hi/ 

18.55 

< 

CO2 (11102) 

-4 

002(10002) 

4 

hi/ 

1.96 

k 

CO2 (11101) 

-4 

002(10002) 

4 

hi/ 

0.117 

k 

002(00011) 

002(10002) 

4 

hi/ 

0.467 

k 

002(10012) 

-4 

002(10002) 

4 

hi/ 

411.14 

k 

17.18 


TABLE  Al.  (CONTINUED) 


Reaction 

Rate  Constant* 

Reference 

C02( 20012) 

C02( 10002) 

hu 

22.28 

k 

C02( 11102) 

002(02201) 

+ 

hu 

0.500 

k 

C02(03301) 

-4 

C02(02201) 

hi/ 

4.58 

k 

CO2 ( lllOl) 

CO2 (02201) 

+ 

hi/ 

1.17 

k 

C02(02211) 

C02(02201) 

hi/ 

414.50 

k 

002(11102) 

C02( 10001) 

+ 

hi/ 

0.028 

k 

CO2 (11101) 

CO2(10001) 

+ 

hi/ 

2.46 

k 

002(00011) 

-> 

002(10001) 

+ 

hi/ 

0.442 

k 

CO2 ( 10011) 

002(10001) 

4 

hi/ 

408.85 

k 

C02( 20012) 

-# 

002(10001) 

+ 

hi/ 

14.99 

k 

002(20011) 

-4 

002(10001) 

+ 

hi/ 

30.06 

k 

002(11112) 

-4 

002(11102) 

4 

hi/ 

404.14 

k 

CO2 (03311) 

-4 

002(03301) 

4 

hi/ 

408.56 

k 

002(11111) 

002(11101) 

4 

hi/ 

400.33 

k 

Rate  constants  are  in  units  of  1/s  for  unimolecular  reactions, 

cm^/molecule-s  for  bitnolecular  reactions  and  cm^/molecule-s  for 
termolecular  reactions. 

^  S.  L.  Lunt,  C.  T.  Wickham- Jones ,  and  C.  J.  S.  M.  Simpson,  "Rate 

Constants  for  the  Deactivation  of  the  15  jum  Band  of  Carbon  Dioxide 
by  the  Collision  Partners  CH3F,  CO2,  N2,  Ar  and  Kr  over  the 
Temperature  Range  300  to  150  K, "  Chem.  Phvs.  Letters.  115,  60 
(1985).  C.  J.  S,  M.  Simpson,  T.  R.  D.  Chandler,  "A  Shock  Tube 
Study  of  Vibrational  Relaxation  in  Pure  CO2  and  Mixtures  of  CO2 
with  the  Inert  Gases,  Nitrogen,  Deuterium  and  Hydrogen,"  Proc.  Rov. 
Soc.  Lond.  A..  317.  265  (1970. 

^  The  value  of  the  deexcitation  rate  constant  given  here  is 

approximately  a  factor  of  three  lower  than  that  reported  by  R.  D. 
Sharma  and  P.  P.  Wintersteiner ,  "Role  of  Carbon  Dioxide  in  Cooling 
Planetary  Thermospheres,"  Geophvs.  Res.  Letters.  17.  2201  (1990). 

R.  E.  Center,  "Vibrational  Relaxation  of  CO2  by  O  Atoms,"  J.  Chem. 
Phvs. .  59,  3523  (1973). 

c  Estimate  based  on  SSH  theory,  J‘nin2ln3X-nin2-ll-ln3X  = 

9nin2-ll-ln3x"2>^01101-*00001- 
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®  Deexcitation  rate  constants  for  Fermi  coupled  states  are  assumed  to 
be  equal. 

®  C.  Dang,  J.  Reid,  and  B.  K.  Garside,  "Dynamics  of  the  CO2  Lower 

Laser  Levels  as  Measured  with  a  Tunable  Diode  Laser,"  AppI.  Phvs. 
B,  31,  163  (1983). 

f  Estimate  based  on  >^ni  +  ln2ln3X-+nin2+21+2n3X '  °  ^1000x->0220x  •  ' 

^n3n2ln3X-*n3n2ln3X  '  “  ^1000x-*1000x  ‘  • 

9  W.  A.  Rosser,  Jr.,  A.  D.  Wood,  and  E.  T.  Gerry,  "Deactivation  of 
Vibrationally  Excited  Carbon  Dioxide  (1^3)  by  Collisions  with  Carbon 
Dioxide  or  with  Nitrogen,"  J.  Chem.  Phvs..  50.  4996  (1969). 

h  J.  H.  W.  Cramp  and  J.  D.  Lambert,  "Vibrational  Relaxation  of 

C02(i'3)  by  O  Atoms,"  Chem.  Phvs.  Letters.  22,  146  (1973). 

i  Estimate  baaed  on  koooii-.ni^n2ln3x’  =  >^00011-03301®"^®^^’’- 

3  Estimate  based  on  k„j^n2ln3X-*n3n2ln3-lx  >‘00011-00001- 

>^  L.  S.  Rothman,  et  al.,  "The  HITRAN  Database;  1986  Edition,"  AppI . 
Opt. .  26,  4058  (1987). 
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TABLE  A2.  REACTIONS  AND  RATE  CONSTANTS  FOR  O3  MODEL. 


Reaction 


Rate  Constant*  Reference 


0  +  O2  +  M  -»  03(000)  + 

M 

2.79xl0~35e510.0/T 

a,  b 

0  +  O2  +  M  -»  03(001)  + 

M 

2.35xl0“35e510.0/T 

a,  b 

0  +  O2  +  M  -*  O3  ( 002 )  + 

M 

1. 85x10-35^510. 0/T 

a,b 

0  +  O2  +  M  -*  03(003)  + 

M 

1.42xl0-35e^^°-°/'^ 

a,  b 

0  +  O2  +  M  -♦  03(004)  + 

M 

1.06xl0-35e^l0-0/T 

a,b 

0  +  O2  +  M  -»  03(005)  + 

M 

7.59xl0-36e510-0/T 

a,  b 

0  +  O2  +  M  -♦  03(006)  + 

M 

4.73xl0-36e^^°-°/’^ 

a,b 

0  +  O2  +  M  -♦  03(007)  + 

M 

2.53xl0-36e^^°*°/’' 

a,b 

0  +  O2  +  M  -»  03(008)  + 

M 

8.80xl0-37e510-0/T 

a,b 

0  +  O2  +  M  -»  03(009)  + 

M 

2.20xl0-37e510-0/T 

a,b 

03(010)  +  M  03(000) 

+  M 

1.73x10-15x1/2 

c 

03(010)  +  0  a  03(000) 

+  0 

1.73x10-13x1/2 

d 

03(020)  +  M  a  03(010) 

+  M 

3.46x10-15x1/2 

e 

03(020)  +  0  a  03(010) 

+  0 

3.46x10-13x1/2 

e 

03(100)  >  M  a  03(010) 

+  M 

2.89x10-15x1/2 

c 

03(100)  +  0  a  03(010) 

+  0 

5.20x10-13x1/2 

d 

03(100)  +  M'  a  03(001) 

+  M' 

1.00x10-11 

b,  c 

03(200)  +  M  a  03(100) 

+  M 

5.80x10-15x1/2 

e 

03(200)  +  0  a  03(100) 

+  0 

1.04x10-12x1/2 

e 

03(200)  +  M'  a  03(101) 

+  M' 

1.00x10-11 

f 

03(001)  +  M  a  03(010) 

+  M 

2.89x10-15x1/2 

c 

03(001)  +  0  a  03(010) 

+  0 

5.20x10-13x1/2 

d 

03(002)  +  M  a  03(001) 

+  M 

5.80x10-15x1/2 

e 

03(002)  +  0  a  03(001) 

+  0 

1.04x10-12x1/2 

e 

03(003)  +  M  a  03(002) 

+  M 

8.70x10-15x1/2 

e 

03(003)  +  0  a  03(002) 

+  0 

1.56x10-12x1/2 

e 

03(004)  M  a  03(003) 

+  M 

1.16x10-14x1/2 

e 

03(004)  +  0  a  03(003) 

+  0 

2.08x10-12x1/2 

e 

03(005)  +  M  a  03(004) 

+  M 

1.45x10-14x1/2 

e 
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TABLE  A2.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


03(005) 

+ 

0 

03(004) 

+ 

0 

03(006) 

+ 

M 

03(005) 

+ 

M 

03(006) 

0 

- 

03(005) 

+ 

0 

03(007) 

+ 

M 

03(006) 

M 

03(007) 

+ 

0 

- 

03(006) 

+ 

0 

03(008) 

+ 

M 

a 

03(007) 

+ 

M 

03(008) 

+ 

0 

- 

03(007) 

+ 

0 

03(009) 

+ 

M 

= 

03(008) 

+ 

M 

03(009) 

+ 

0 

03(008) 

+ 

0 

03(001) 

+ 

0 

-» 

O2  t  O2 

03(002) 

+ 

0 

-» 

O2  +  O2 

03(003) 

+ 

0 

-* 

O2  +  O2 

03(004) 

+ 

0 

-» 

O2  +  O2 

03(005) 

+ 

0 

-» 

O2  +  O2 

03(006) 

+ 

0 

-* 

O2  +  O2 

03(007) 

+ 

0 

-♦ 

O2  +  O2 

03(008) 

+ 

0 

-» 

O2  +  O2 

03(009) 

+ 

0 

-* 

O2  +  O2 

03(011) 

+ 

M 

S! 

03(010) 

+ 

M 

03(011) 

+ 

0 

03(010) 

+ 

0 

03(011) 

+ 

M 

- 

03(001) 

+ 

M 

03(011) 

+ 

0 

03(001) 

+ 

0 

03(110) 

+ 

M 

=: 

03(010) 

+ 

M 

03(110) 

+ 

0 

03(010) 

0 

03(110) 

+ 

M 

= 

03(100) 

+ 

M 

03(110) 

+ 

0 

s 

03(100) 

+ 

0 

03(101) 

+ 

M 

s 

03(001) 

+ 

M 

03(101) 

+ 

0 

= 

03(001) 

+ 

0 

03(101) 

+ 

M 

03(100) 

+ 

M 

2.60x10"12t1/2 

e 

1.74x10“1'*t1/2 

e 

3. 12x10“I2t1/2 

e 

2.03x10"14t1/2 

e 

3.64x10"12t1/2 

e 

2.32x10"1‘*t1/2 

e 

4.16x10-12^1/2 

e 

2.61x10“14t1/2 

e 

4.68x10-12x1/2 

e 

1. 00x10-11 

b,  d 

1.00x10-11 

b,d 

1. 00x10-11 

b,  d 

1.00x10-11 

b,  d 

1.00x10-11 

b,  d 

1.00x10-11 

b,  d 

1.00x10-11 

b,d 

1.00x10-11 

b,  d 

1.00x10-11 

b,d 

2.89x10-15x1/2 

g 

5.20x10-13x1/2 

g 

1.73x10-15x1/2 

g 

1.73x10-13x1/2 

g 

2.89x10-15x1/2 

g 

5.20x10-13x1/2 

g 

1.73x10-15x1/2 

g 

1.73x10-13x1/2 

g 

2.89x10-15x1/2 

g 

5.20x10-13x1/2 

g 

2.89x10-15x1/2 

g 
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TABLE  A2.  (CONTINUED) 


Reaction 

Rate  Constant* 

Reference 

03(101) 

+  O  =  03(100)  + 

0 

5.20x10-13t1/2 

g 

03(101) 

+  M'  =  03(002) 

*■  M* 

1.00x10-11 

f 

03(111) 

+  M  =  03(101)  + 

M 

1.73x10-15t1/2 

g 

03(111) 

+  O  =  03(101)  + 

0 

1.73x10-12x1/2 

g 

03(111) 

+  M  =  03(011)  + 

M 

2.89x10-15x1/2 

g 

03(111) 

+  O  =  03(011)  + 

0 

5.20x10-13x1/2 

g 

03(111) 

+  M  =  03(110) 

M 

2.89x10-15x1/2 

g 

03(111) 

+  0  =  03(110)  + 

0 

5.20x10-13x1/2 

g 

03(010) 

-*  03(000)  +  hi/ 

0.224 

h 

03(001) 

-»  03(000)  +  hi/ 

11.00 

h 

03(100) 

-♦  03(000)  +  hi/ 

0.594 

h 

03(011) 

-*  03(000)  +  hi/ 

0.116 

h 

03(110) 

03(000)  +  hi/ 

0.053 

h 

03(002) 

-♦  03(000)  +  hi/ 

0.342 

h 

03(101) 

-»  03(000)  +  hi/ 

3.67 

h 

03(200) 

-*  03(000)  +  hi/ 

0.106 

h 

03(111) 

-»  03(000)  +  hi/ 

0.131 

h 

03(003) 

-*  ©3(000)  hi/ 

0.746 

h 

03(020) 

-*  ©3(010)  +  hi/ 

0.446 

h 

03(011) 

-*  ©3(010)  +  hi/ 

10.39 

h 

03(110) 

-•  ©3(010)  +  hi/ 

0.291 

h 

03(111) 

-*  ©3(010)  +  hi/ 

3.67 

h 

03(.  )2) 

-*  ©3(001)  +  hi/ 

20.72 

h 

03(101) 

-♦  ©3(001)  +  hi/ 

0.437 

h 

03(111) 

-♦  ©3(001)  +  hi/ 

0.053 

h 

03(003) 

->  ©3(001)  +  hi/ 

1.03 

h 

03(101) 

-♦  ©3(100)  +  hi/ 

9.68 

h 

03(200) 

-*  ©3(100)  +  hi/ 

1.18 

h 

03(111) 

-*  O3  ( 100 )  +  hi/ 

0.116 

h 
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TABLE  A2.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


03(111) 

-♦  03(011) 

+ 

hi/ 

0.593 

h 

03(111) 

-+  ©3(110) 

+ 

hi/ 

11.01 

h 

03(003) 

-»  ©3(002) 

+ 

hi/ 

30.44 

h 

03(004) 

-►  ©3(003) 

+ 

hi/ 

36.34 

h 

©3(005) 

-»  ©3(004) 

+ 

hi/ 

43.08 

h 

©3(006) 

-*  ©3(005) 

+ 

hy 

48.49 

h 

03(007) 

■>  ©3(006) 

+ 

hy 

52.71 

h 

03(008) 

-»  O3  ( 007  ) 

+ 

hy 

56.12 

h 

03(009) 

-♦  ©3(008) 

+ 

hy 

58.67 

h 

Rate  constants  are  in  units  of  1/s  for  unimolecular  reactions, 
cm^/molecule-s  for  bimolecular  reactions  and  cm^/molecule-s  for 
termolecular  reactions. 

^  R.E.  Huie,  J.T.  Herron,  and  D.D.  Davis,  "Absolute  Rate  Constants 
for  the  Reaction  of  0+02+M-K)3+M  over  the  Temperature  Range  200-346 
K,"  J.  Phvs.  Chem..  76,  2653  (1972). 

^  W.T.  Rawlins,  "Chemistry  of  Vibrationally  Excited  Ozone  in  the 
Upper  Atmosphere,"  J.  Geophvs.  Res..  90.  12283  (1985). 

^  S.M.  Adler-Golden  and  J.I.  Steinfeld,  "Vibrational  Energy  Transfer 
in  Ozone  by  Infrartd-Ultraviolet  Double  Resonance,"  Chem.  Phvs. 
Letters .  76.  479  (1980). 

G.A.  West,  R.E.  Weston,  Jr.,  and  G.W.  Flynn,  "Deactivation  of 
Vibrationally  Excited  ozone  by  0(3p)  Atoms,"  Chem.  Phvs.  Lett..  42 , 
488  (1976). 

®  Estimate  based  on  SSH  theory,  k^-jn-i  =  nk^-^Q. 
f  Estimate  based  on  kiinn-,i_in,n+i  =  kioo-»001- 

9  Estimate  based  on  =  kQjo_,QQQ,  =  k3Q0-*000' 

>tlmn-»lmn-l  =  kooi->ooo* 

^  L.  S.  Rothman,  et  al.,  "The  HITRAN  Database:  1986  Edition,"  Appl . 
Qpt. .  26f  4058  (1987) . 
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TABLE  A3 


REACTIONS  AND  RATE  CONSTANTS  FOR  NO  MODEL 


Reaction 

Rate  Constant* 

Reference 

N0(1) 

+ 

O2  ^  NO(0)  -*•  O2 

2.6x10-14 

a 

NO(2) 

■f 

O2  ^  NO(l)  +  O2 

6.6x10-14 

a 

NO(l) 

+ 

0  =  NO(0)  +0 

6.5x10-11 

b 

NO(2  ) 

+ 

0  NO(l)  +0 

2.17x10-11 

c 

NO(2) 

+ 

0  =  NO(0)  +  0 

4.33x10-11 

c 

NO(l) 

NO(0)  +  hi/ 

13.30 

d 

NO(2  ) 

NO(l)  +  hu 

25.43 

d 

NO(2) 

-* 

NO(0)  +  hi/ 

0.857 

d 

Rate  conatants  are  in  units  of  cm^/tnolecule-sec  for  biroolecular 
processes  and  1/sec  for  the  Einstein  A  coefficients. 


B.  D.  Green,  B.  E.  Caledonia,  R.  E.  Murphy,  and  F.  X.  Robert,  "The 
Vibrational  Relaxation  of  NO(v=l-7)  by  O2,”  J.  Chem.  Phvs..  76, 
2441  (1982). 


R.  P.  Fernando  and  I.  W.  M.  Smith,  "Vibrational  Relaxation  of  NO  by 
Atomic  Oxygen,"  Chem.  Phvs.  Letters.  66.  218  (1979). 


K.  Glanzer  and  J.  Troe,  "Vibrational  Relaxation  of  NO  in  Collisions 
with  Atomic  Oxygen  and  Chlorine,"  J.  Chem.  Phvs..  63 .  4352  (1975). 
For  an  estimate  of  the  importance  of  multiguantum  quenching,  see  M. 
Quack  and  J.  Troe,  "Complex  Formation  in  Reactive  and  Inelastic 
Scattering;  Statistical  Adiabatic  Channel  Model  of  Unimolecular 
Processes  III,"  Ber.  Bunsenoes.  Phvsik  Chem.,  79 .  170  (1975). 


d 


L.  S.  Rothman,  et  al.,  "The  HITRAN  Database: 
Opt. .  26,  4058  (1987). 


1986  Edition,"  Appl . 
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TABLE  A4.  REACTIONS  AND  RATE  CONSTANTS  FOR  H2O  MODEL. 


Reaction 

Rate  Constant* 

Reference 

H20(010) 

+ 

M  =  H20(000)  +  M 

1.10xl0-10e-54.0/Tl/3 

a 

H20(100) 

+ 

N2  «  H20(020)  +  N2 

2.7x10~14t1/2 

c 

H20(100) 

O2  -  H20(020)  +  O2 

1.9x10'14t1/2 

c 

H20(100) 

+ 

0  ^  H20(020)  +  0 

4x10-12 

d 

H20(001) 

+ 

N2  =  H20(020)  +  N2 

2.7x10“14t1/2 

c 

H20(001) 

+ 

O2  =  H20(020)  +  O2 

1.9x10-14t1/2 

c 

H20(001) 

+ 

0  =  H20(020)  +  0 

4x10-12 

d 

H20(001) 

+ 

M  -  H20(100)  +  M 

1.4x10-11t1/2 

e 

H20(020) 

+ 

M  s  H20(010)  +  M 

2.20xlO-10e-^4.0/Tl/3 

f 

H20(030) 

+ 

M  «  H20(020)  +  M 

3.30xlO-10e-54.0/Tl/3 

f 

H20(110) 

+ 

M  «  H20(100)  +  M 

1.10xlO-10e-54.0/Tl/3 

a 

H20{011) 

+ 

M  «  H20(001)  +  M 

1. 10xl0-10e”^'^'°'^1’^^^ 

a 

H20(110) 

+ 

N2  «  H20(030)  +  N2 

2.7x10-14t1/2 

g 

H20(110) 

O2  =  H20(030)  +  O2 

1.9x10-14t1/2 

g 

H20(011) 

+ 

N2  »  H20(030)  +  N2 

2.7x10-14x1/2 

g 

H20(011) 

+ 

O2  =  H20(030)  +  O2 

1.9x10-14x1/2 

g 

H20(010) 

02(0)  =  H20(000)  +  02(1) 

2.2x10-11/x1/2 

a 

H20(020) 

02(0)  =  H20(010)  +  02(1) 

2.2x10-11/t1/2 

a 

H20(030) 

+ 

02(0)  «  H20(020)  +  02(1) 

2.2x10-11/x1/2 

a 

H20(110) 

+ 

02(0)  =  H20(100)  +  02(1) 

2.2x10-11/x1/2 

a 

H20{011) 

+ 

02(0)  =  H20(001)  +  02(1) 

2.2x10-11/x1/2 

a 

H20(010) 

H20(000)  +  hi/ 

20.46 

h 

H20(020) 

H20(000)  +  hi/ 

0.501 

h 

H20(100) 

“4 

H20(000)  +  hi/ 

5.79 

h 

H20(001) 

-♦ 

H20(000)  +  hi/ 

82 . 51 

h 

H20(030) 

-4 

H20(000)  +  hi/ 

0.005 

h 

H20(110) 

-4 

H20(000)  +  hi/ 

0.596 

h 

H20(011) 

-4 

H20(000)  +  hi/ 

17.09 

h 

H20(020) 

H20(010)  +  hi/ 

35.27 

h 
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TABLF  A4.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


H20(100)  -*  H20(010)  +  hi/  1.34 

H20(001)  -♦  H20(010)  +  hu  2.21 

H20(030)  -♦  H20(010)  +  hi/  1.13 

H20(110)  -»  H20(010)  +  hi/  5.12 

H20(011)  -»  H20(010)  +  hi/  77.20 

H20(030)  -*  H20(C20)  +  hi/  40.28 


Rate  constants  are  in  units  of  1/s  for  unimolecular  reactions, 
ctn^/molecule-s  for  bimolecular  reactions  and  cm^/molecule-s  for 
termolecular  reactions. 

H.  E.  Bass  and  F.  D.  Shields,  "Vibrational  Relaxation  and  Sound 
Absorption  in  O2/H2O  Mixtures,"  J.  Acoust.  Soc.  Am..  56.  856 

(1974) . 

R.  T.  V.  Rung  and  R.  E.  Center,  "High  Temperature  Vibrational 
Relaxation  of  H2O  by  H2O,  He,  Ar,  and  N2,”  J.  Chem.  Phvs..  62 .  2187 
(1975) . 

J.  Finzi,  F.  E.  Hovis,  V.  N.  Panfilov,  P.  Hess,  and  C.  8.  Moore, 
"Vibrational  Relaxation  of  Water  Vapor,"  J.  Chem.  Phvs..  67 .  4053 
(1977) . 

P.  F.  Zittel  and  D.  E.  Masturzo,  vibrational  Relaxation  of  H2O  from 
295  to  1020  K,"  J.  Chem.  Phvs..  90,  977  (1989). 

Gas  kinetic  estimate. 

Estimate  based  on  SSH  theory,  =  nkj-^g* 

Estimate  based  on  kin0-»0n+20  =  '^100-*020  *^0nl-*0n+20  =  •'001-»020- 

L.  S.  Rothman,  et  al.,  "The  HITRAN  Database:  1986  Edition,"  AppI . 

Opt. .  26.  4058  (1987). 


TABLE  A5.  REACTIONS  AND  RATE  CONSTANTS  FOR  CO  MODEL 


Reaction 

Rate  Constant* 

Reference 

CO(l) 

+ 

M  =  CO(0)  +  M 

9.90xl0“9e“^^®- 

a 

CO(2) 

+ 

M  =  CO(l)  +  M 

1.98xl0"8e"l®®- 

b 

CO(l) 

+ 

O  e  CO(0)  +  O 

2.82xl0-9e-75.4/Tl/3 

c,d 

CO(2) 

+ 

0  =  CO(l)  +  O 

2.82xl0-9e-75.4/Tl/3 

e 

CO(2) 

+ 

O  «  CO(0)  +  O 

2.82xl0-9e-^5.4/Tl/3 

e 

CO(0) 

+ 

N2(1)  -  CO(l)  +  N2(0) 

6.98xl0“13e-25.6/T3/3 

a 

CO(l) 

C0(0)  +  hu 

30.96 

f 

CO(2) 

CO(l)  +  hf 

60.45 

f 

CO(2) 

CO(0)  +  hi/ 

1.03 

f 

Rate  constants  are  in  units  of  1/s  for  unimolecular  reactions, 
cm^/molecule-s  for  bimolecular  reactions  and  cm^/molecule-s  for 
termolecular  reactions. 

R.L.  Taylor,  "Energy  Transfer  Processes  in  the  Stratosphere,"  Can. 
J.  Chem..  1436  (1974). 

Estimate  based  on  SSH  theory,  k2i  = 

R.  E.  Center,  "Vibrational  Relaxation  of  CO  by  O  Atoms,"  J.  Chem. 
Phvs. .  58f  5230  (1973). 

M.  E.  Lewittes,  C.  C.  Davis,  R.  A.  Mcfarlane,  "Vibrational 
Deactivation  of  CO  by  Oxygen  Atoms,"  J.  Chem.  Phvs..  69,  5230 

(1978) . 

Estimate. 

L.  S.  Rothman,  et  al.,  "The  HITRAN  Database:  1986  Edition,"  AppI . 
Opt. .  26,  4058  (1987). 
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